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ABSTRACT 


This document provides a general survey of the field of space-vehicle structural 
vibration, which is induced by acoustic and aerodynamic noise and certain 
mechanical excitation. The fundamental sources and mechanisms of the vibration 
are identified; analytical methods of vibration prediction are described. 
Laboratory, field, and flight testing are discussed, including problems of 
qualification- and acceptance- test selection. Although the advantages and 
limitations of the various analytical and experimental methods are assessed, 
specific methods are seldom recommended because selection of the methods is 
often influenced by the circumstances of the particular application, and by 
economic and scheduling factors. 

A list of program objectives is furnished outlining a general procedure for 
ensuring the structural integrity and operational performance of launch vehicles 
and spacecraft exposed to structural vibration. Then, to ensure that these 
objectives are satisfied, general recommendations are made for implementing the 
various design, analysis, and testing techniques during various phases of 
vehicle development. 
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ASSESSMENT OF SPACE VEHICLE AEROACOUSTIC- VIBRATION 


PREDICTION, DESIGN, AND TESTING 

By Harry Himelblau, Space Division of North American Rockwell 

Corporation; C. Myron Fuller, McDonnell Douglas Astronautics 

Company - Western Division; and Terry D. Scharton, Bolt 

Beranek and Newman Inc. 

1 . INTRODUCTION 

Various load and environmental conditions are applied to a space vehicle, its 
subsystems, or assemblies during a mission - conditions that must be considered 
in design and test to avoid mission failure or degradation. Mechanical vibration 
of the structure is an extremely important condition to be considered, since 
it may cause overstress or fatigue of the structure and equipment, or may 
cause malfunctions of the equipment. This report restricts itself to consider- 
ation of the integrity of structure; the performance and structural integrity 
of equipment will not be discussed herein. This document does however include 
a discussion of vibration design and test requirements for this equipment. 

Structural vibration is often most predominant during the launch -and- as cent 
phase of flight, and, depending upon the mission parameters and the structural 
configuration, may be critical during space flight and/or atmospheric entry 
and flight. For example, figure 1 shows the time history of the instantaneous 
vibration and the time-averaged rms vibration of a point on the S-IVB stage of 
the uprated Saturn I space vehicle, measured during the launch -and -as cent 
phase. Significant vibration is observed during liftoff (T + 0 sec), the 
transonic period (T + 47) , the supersonic period near the occurrence of the 
maximum aerodynamic pressure (1 + 63) , and during S-IVB powered flight 

(from T + 147 to T + 434). Vibration spectra, in the form of acceleration 
spectral density vs frequency plots, are also shown for liftoff, transonic, 
and two periods during the S-IVB powered flight. Depending upon the structural 
and aerodynamic configuration of the vehicle and the location of the vibration- 
monitoring points on the structure, the relative magnitude of the vibration 
during these periods will vary considerably from that shown in figure 1 . Loca- 
tions near rocket engines, for example, may undergo high vibration during 
propulsion operation, and relatively low vibration at other time periods. 
Locations near the forward end of the vehicle, such as spacecraft structures, 
may undergo high vibration during transonic and/or supersonic periods, and 
relatively low vibrations at other times. Before flight, vibration may be 
applied to the vehicle, its subsystems, or assemblies during transportation 
from the manufacturer to the launch pad, or during the static firing of propul- 
sion subsystems. Table I lists the mission phases and time periods where 
significant vibration may possibly be observed. 

Usually, the attempt is made to integrate these environments into the vehicle 
development program by means of design and test requirements. The selection 
of these requirements can be vital to mission success and cost, and is there- 
fore one of the subjects of this report. Requirements which are too low. 
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TABLE I. -SIGNIFICANT STATIONARY VIBRATION ENVIRONMENTS 


Mission phase 

Source of vibration 

Type, source, and magnitude 
of vibration 

Prelaunch 

Acceptance testing 

Random and/or periodic vibration 
from laboratory vibration and/or 
acoustic test, if used. Often 
high magnitude. 


Transportation 



Land 

Random vibration from truck or 
inplant dolly. Periodic vibration 
also possible. Random and 
periodic vibration from railroad 
transportation. Usually (but not 
always) low magnitude and long 
duration. 


Air 

Random vibration from jet aircraft. 
Periodic vibration also possible 
adjacent to jet engines. Periodic 
vibration from propeller aircraft. 
Random and periodic vibration from 
turboprop aircraft. Random 
vibration from gust loading on 
all aircraft at low frequencies. 
Usually (but not always) low 
magnitude and long duration. 


Sea 1 

Random and periodic vibration 
from ships. Almost always low 
magnitude (except near fantail) 
and long duration. 


Static firing 

Random vibration from engine - 
induced acoustic noise. Often 
high magnitude. Directly trans- 
mitted random and/or periodic 
vibration from engines, turbo- 
pumps , and auxiliary equipment . 
Often high magnitude near source. 


Ground wind loads 

Random and/or periodic vibration 
at low frequencies. Large limit - 
amplitude self-excited oscillation 
may exceed limit stresses if 
system is unstable. 




TABLE I. -SIGNIFICANT STATIONARY VIBRATION ENVIRONMENTS - Continued 


Mission phase Source of vibration 


Launch and 
ascent 


Liftoff 

Wind loads - 
engine interaction 

Aeroelastic 
interaction 
and flutter 


Pogo interaction 


Propellant sloshing 


Transonic buffeting 


Supersonic flight 


Booster or upper- 
stage operation 


Type, source and magnitude 
of vibration 

Same as 4 ‘static firing,* * except 
vibration spectra may be different. 

Transient or periodic vibration at 
low frequencies. Divergent 
vibration if system is unstable. 

Transient or periodic vibration 
of structural assemblies exposed 
to aerodynamic flow. Divergent 
motion or limit amplitude if 
system is unstable. 

Transient or periodic vibration 
along the longitudinal vehicle 
axis, from engine- thrust oscilla- 
tions, causing structural vibration 
of propulsion assemblies, and in 
turn causing propellant-flow and 
engine- thrust oscillations. 
Divergent motion or limit ampli- 
tude if system is unstable. 

Transient or periodic vibration 
at low frequencies. Divergent if 
slosh-structure-control system is 
unstable. 

Random vibration from aerodynamic 
turbulence, sometimes c ‘ amplified * * 
by unstable shock waves. May 
cause local or overall vehicle 
vibration. Often high magnitude, 
usually at forward end of vehicle, 
and sometimes transient. 

Random vibration from aerodynamic 
turbulence in relatively stable 
q^ ax boundary layer. Often high 

magnitude, usually at forward end 
of vehicle. 

Directly transmitted random and/or 
periodic vibration from engines, 
turbopumps , and auxiliary 
equipment. Often high magnitude 
near source. 


TABLE I. -SIGNIFICANT STATIONARY VIBRATION ENVIRONMENTS - Concluded 


Mission Phase 

Source of vibration 

Type, source and magnitude 
of vibration 

Launch and 

Abort 

If abort rocket is aft, same as 

ascent 


‘ Booster or upper-stage opera- 

(concluded) 


tion.’’ If abort rocket is 
forward, random vibration from 
exhaust flow. Possible inter- 
action with or increased 
transonic buffeting or boundary 
layer turbulence. Often high 



magnitude. 

Space flight. 

Booster or upper- 

Same as under * ‘Launch and 


stage operation 

ascent . * 5 

Entry and 

Hypersonic flight 

Random vibration from aerodynamic 

atmospheric 


turbulence. Usually lower 

flight 


magnitude due to low density of 
boundary layer and, in some cases, 
separated base flow. 


Deceleration loads 

Random vibration from retro - 
rocket exhaust flow. Possible 
interaction with hypersonic 
turbulence. Random vibration 
from aerodynamic turbulence on 
parachute. Usually low magnitude. 


compared to the vehicle environment, may jeopardize the mission because of 
failure to discover design or fabrication inadequacies. On the other hand, 
excessive requirements may obtain unnecessarily high reliability, and will 
usually cause excessive costs, schedule slippage, and/or excessive vehicle 
weight. 

By definition, structural vibration is an oscillatory motion of a mechanical 
system relative to a frame of reference, and is described by certain param- 
eters which are considered to be stationary or steady state. This vibration 
may be periodic or random, or both. In space vehicles, this vibration occurs 
only for definable periods of time, and therefore it is theoretically tran- 
sient. However, there are several phases of the mission where the significant 
parameters defining the vibration (such as the rms value) are invariant or are 
slowly changing with time, as observed in figure 1. During these periods, the 
effects of the vibration on vehicle integrity and performance may be determined 
adequately through the use of stationary or steady-state assumptions in the 


vibration analysis and test. The minimum period that a vibration may be 
considered as stationary is usually dependent on the characteristics of the 
source, the structural characteristics (especially the damping), the failure 
and malfunction parameters, the accuracy desired, and sometimes common aerospace 
practice or the investigator’s judgment. Vibration not considered sufficiently 
stationary or steady state should be considered as transient, and is not dis- 
cussed in this report. 

The sources of space-vehicle vibration discussed in this report are acoustic 
noise, aerodynamic noise, and mechanically induced excitation. The predominant 
frequencies of these sources usually cover a wide range, often from 1 0 Hz to 
10 kHz. In most locations on the vehicle, however, vibration above 2 kHz 
is seldom a cause of failure, so that the main emphasis is devoted to the 
range from 10 Hz to 2 kHz. In this frequency range, subsections of the vehicle 
often vibrate independently of other subsections (e.g., entire stages, individual 
skin panels, longerons, stiffeners, or frames, or equipment items on structure). 
At lower frequencies, the vehicle usually vibrates 4 ‘as a whole’’ in its various 
lateral, longitudinal, and torsional modes. Often, this behavior begins to 
disappear before the 10th lateral mode, whose resonant frequency usually varies 
inversely with the size of the vehicle. In the larger space vehicles, this 
resonant frequency usually occurs below 25 Hz. 

It should be emphasized that the initial design of the space vehicle is almost 
never made considering structural vibration. Most vehicle structures are 
initially designed to static and sometimes thermal loads, with some overdesign 
factors specified for such additional loads as structural vibration. It will 
be necessary to ascertain whether enough overdesign has been incorporated to 
obviate the need to redesign, considering structural vibration and other loads 
and environments. 
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SYMBOLS 


A 


B p (x,x) 

[C] 

C p (5,S',f) 

c r (f) 


c 


C L 


D 

Q 

E 

F, F(x,t) 


f 

f 


f., f 
i n 




area of structure exposed to spatially distributed 

2 

applied loading, in. 
radius, in. 

transfer function between internal force or 
moment at location x and applied force at location x 

damping matrix, lb -sec/in. 

cospectral density function of fluctuating pressures at 

2 

frequency f and locations £ and £ ! , (psi) /Hz 

ratio of vibration response to applied loading at frequency 
f 

specific heat (subscripts denote medium), in.^/sec^ °F or 
Btu in. /lb sec^ °F 

speed of sound in acoustic medium, in. /sec 

speed of longitudinal waves in a plate or shallow shell 
(subscripts denote structural system), in. /sec 

fatigue damage; vehicle diameter, in.; flexural rigidity, 
in. -lb 

inplane rigidity, lb/in. 

Young’s modulus, psi 

instantaneous external force (subscripts denote single 
amplitude or source) , lb 

frequency, Hz 

acoustic-structure coincidence frequency, Hz 

th 

natural frequency in i mode, Hz 
ring frequency, Hz 

cross-spectral density functions for sources I and K at 
frequency f 
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G p (5,?',f) 


G (x,f) 

G (X,x' ,f) 
g 

H.(f) 


h(x) 


I 

sp 

I(u,q) 


J 


2 

ik 


(f) 


[K] 

K d 

K f 

% 

CK 


k 


cross-spectral density function of fluctuating pressures 

at frequency f and locations £ and £’ of the structure 

2 

exposed to the pressure field, (psi) /Hz 

reference pressure spectral density function at frequency 
f, (psi) 2 /Hz 

spectral density function of parameter at location x and 

2 

frequency f (subscript denotes parameter), (...) /Hz 

cross -spectral density function of parameter at locations x 

2 

and x ! (subscripts denote parameter), (...) /Hz 

2 

acceleration due to earth gravity, in ./sec 

th 

frequency response function in i mode (asterisk denotes 
complex conjugate) 

transfer function between any response parameter at location 
x and any applied loading parameter I (asterisk denotes 
complex conjugate) 

thickness of structure at location x (subscripts denote core 
half- thickness and face-sheet thickness for sandwich shell) , 
in. 

specific impulse of engine, sec 

ratio of incomplete to complete gamma function 

cross -joint acceptance function of pressure field with the 
i^ and k^ structural mode shapes. 

stiffness matrix, lb/in. 

3 

ratio of stress to displacement, lb sec/in. 
stress-concentration factor in fatigue 
mode-shape factor 

sine- to-random fatigue conversion factor 

thermal conductivity (subscripts denote medium), lb/sec °F 
or Btu/in. sec °F 
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% 

M 

[M] 

m 

N 


N. 


n. 

x 


length, in. 

moment -per -unit length (subscripts denote directions), lb 

2 

mass or inertia matrix, lb-sec /in. 

th 2 

modal or generalized mass in i mode, lb sec /in. 

2 

mass of structure, lb-sec */in. 

number of modes; number of measurements; number of engines; 
inplane-f orce-per-unit length (subscripts denote direc- 
tions) , lb/in. 

number of cycles to failure at stress amplitude s_^ 

number of applied cycles at stress amplitude s^; scale 
factor 


P 

o 


p T (3 > 3 ) 

1 — o 


P» p(5,t) 


w 

{Q} 


Q. Q ± . 


n 


Q p (e,c f ,f) 

q 

{q} 




static or absolute ambient pressure, psi 

probability that stress peak s . will exceed threshold 
* 1 
stress s^ during exposure time T 

instantaneous pressure at location £, psi 

probability density of stress peaks, (psi) ^ 

t 

column matrix of applied generalized forces in i normal 
modes, lb or in. lb 

th 

quality factor or resonant magnification in i mode; 
transverse f orce-per-unit length (subscripts denote 
directions), lb/in. 

quad-spectral density function of fluctuating pressures 

2 

at frequency f and locations £ and £*, (psi) /Hz 

2 

aerodynamic pressure, lb/ft 

th 

column matrix of coordinate displacements in i normal 
modes, in. or radians 

cross- correlation function of fluctuating pressures at 

2 

locations £ and £ 1 and time delay t, (psi) 
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R 

S 

s 

T 

T r 


rad 


U, u(x,t) 
v, v(x,t) 


W 

W, w (x 5 t ) 


X. 

1 

x,x' jXjX 1 


Z(f) 


3 

r(q 


1 ) 


Y 




radiation resistance of structure in acoustic medium, 
lb-sec/ in. 

2 

total area of structure, in. 
stress (subscripts denote type) , psi 

time at liftoff; exposure time, sec; engine thrust, lb 
time-to-f ailure , sec 
time, sec 

instantaneous displacement in direction at location x 
(subscripts denote partial spatial derivatives), in. 

instantaneous displacement in X^ direction at location x 
(subscripts denote partial spatial derivatives), in. 

weight, lb 

instantaneous displacement in Z direction at location x 
(subscripts denote single amplitude or partial spatial 
derivatives), in. 

rms value of i ^ measurement, g 

locations on structure (underline denotes dummy variable to 
be integrated) 

mechanical driving point impedance of system at frequency f 
distance from neutral plane in Z direction, in. 
logarithmic slope of s^ vs 1SL curve 

coefficient of skewness for statistical distribution 
complete gamma function of parameter q 
mass-attenuation factor 
multimodal-response factor 
viscous damping ratio in i t ^ L mode 
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n 


D 


y 

V 


5’ 


t r . 

i 


P 


a 


P> p(x) 


a 


<f> ± 00 

v» ± (x) 

4^00 


<• 


• > 


SUBSCRIPTS 

a 

cr 

e 


structural damping- loss factor or coefficient; coupling- 
loss factor between systems (subscripts denote structural 
or acoustic systems) 

phase, rad; time-average modal energy for structural and 
acoustic modes, in. -lb 

2 

coefficient of viscosity, lb-sec/in. 

Poission’s ratio (subscripts denote sandwich shell) 

locations on part of structure exposed to spatially 
distributed applied loading (underline denotes dummy 
variable to be integrated) 

dimensionless group of parameters 

2 4 

mass density of acoustic medium, lb-sec /in. 

2 4 

mass density of structure at location x, lb-sec /in. 

standard deviation of parameter (subscript denotes para- 
meter); rms value when mean value is zero or can be 
ignored 

time delay for cross-correlation analysis, sec 

th 

mode shape at location x in Z direction and i mode 

mode shape at location x in direction and i* 1 * 1 mode 

mode shape at location x in direction and i mode 

absolute value of parameter 

spatial average of parameter 

time average of parameter 

dummy variable to be integrated 


acoustic medium (usually air) ; accelerometer 
criteria 

endurance limit; external structure; equipment 
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F 

f 

i 

k 

in 

max 

min 

n 

o 

P 

P 

r 

s 

u 

v 

w 

1 , 2,3 


external force 

fluid, failure 
tti th 

i mode; i measurement; internal structure; peak value 

k *** 1 mode; number of frequency band 

model or subscale structure 

maximum 

minimum 

n *'* 1 mode; new vehicle 
single amplitude 
internal force or moment 

pressure; prototype or full-scale structure 

reference vehicle 

structure 

ultimate; displacement in X^ direction 
displacement in direction 
displacement in Z direction 

system numbers; directions; partial spatial derivatives 


12 



3, FUNDAMENTAL SOURCES OF VIBRATION 


3.1 ACOUSTIC NOISE FROM PROPULS ION -SYSTEM OPERATION 

The fundamental sources of structural vibration are many and varied. Rocket 
engines emit high-velocity exhaust gases which mix with the ambient air, causing 
turbulence in the process. The pressure fluctuations of the turbulence are 
transmitted to the surroundings, including the space vehicle, as acoustic noise 
at the speed of sound. The distributed acoustic field progresses over the 
vehicle surface and causes the structure to vibrate. Since the turbulent mixing 
is a random process, the acoustic noise and the structural vibration are also 
random. Observations show that most of the higher-frequency noise is generated 
by smaller- scale mixing near the engine, while most of the lower- frequency 
noise is generated by larger-scale mixing further downstream. 

If the launch pad or the test stand has a flame deflector, the exhaust gases are 
caused to mix differently with the ambient air than they would if no flame 
deflector were present, and the directivity of the acoustic noise changes 
relative to the surroundings, including the space vehicle, as shown in figure 2. 
Thus, the acoustic field changes, causing a change in the structural vibration. 
The launch pad or test stand, the surrounding terrain, and climatic conditions 
may have similar influences on the acoustic field by providing varying degrees 
of sound reflection and shielding. 



Figure 2 -Rocket-engine acoustic noise sources and directivity patterns to the space vehicle during 
liftoff or static firing. 
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Some space vehicles use multiengine arrangements. These have effects on the 
turbulence generation and the acoustic noise transmission that are different 
from single-engine effects. Rocket-engine noise characteristics are summarized 
in references 1 through 7, which also reference other pertinent documents. 

3.2 MECHANICAL VIBRATIONS FROM PROPULSION- SYSTEM OPERATION 

The mixing of propellant gases inside the rocket-engine combustion chamber 
generates fluctuating pressures (as well as static pressures) on the engine 
wall. The wall vibrations may be transmitted throughout the vehicle, but are 
often rapidly attenuated with distance from the source. In addition, acoustic 
resonances may occur in the combustion chamber. For a liquid-propelled rocket 
engine, these resonances may occur at fixed frequencies, while for a solid- 
propellant rocket motor, these resonances may occur at progressively lower 
frequencies as the propellant burns and the size of the combustion chamber 
increases. These resonances are substantially attenuated or eliminated during 
the design or development of the rocket engine. Failure to do so is one cause 
of engine explosion (refs. 8 through 11). 

Turbopumps and other equipment often have rotating parts operating at high 
speed. Turbine-blade resonance and unbalance of rotating parts may cause 
periodic vibration throughout the vehicle, but this vibration is often rapidly 
attenuated with the increase of distance from the source. 

3.3 AERODYNAMIC FLOW FIELD 

As a space vehicle moves at high speed through the atmosphere, turbulence is 
generated by mixing in the boundary layer. The resulting pressure fluctuation 
is applied over the vehicle surface and causes the structure to vibrate. This 
pressure field, often called aerodynamic noise or buffeting, is transported 
down the vehicle surface at a speed which is somewhat less than the free-stream 
velocity of the vehicle. The speed of the aerodynamic flow seldom equals the 
speed of sound in air. As the vehicle cross-sectional area changes from fore 
to aft, shock waves occur in the aerodynamic flow, causing increases in the 
turbulence and fluctuating pressures. Examples of these shock waves and typical 
effects on the pressure field are shown in figure 3. Note that separated flow 
may occur on a conical section of the vehicle aft of a cylindrical section, and 
vice versa. 

In the transonic regime, the shock waves, which have just formed at various 
locations, are often unstable and tend to move aft. During the launch -and -as cent 
phase, the vehicle speed is increasing in the transonic regime, so that a non- 
stationary or transient pressure field exists, causing a vibration of the 
structure. 

In the supersonic regime, the shock waves are generally stable and less 
susceptible to speed or attitude changes. Thus, the fluctuating pressure field 
and structural vibration are more nearly stationary. Aerodynamic noise 
characteristics are summarized in references 6, 7, and 12 through 20, which 
also reference other pertinent documents. However, much research on aerodynamic 
noise is still needed before this environment can be adequately predicted. 
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SUBSONIC (M =0.90) 

figure 3-Schlieren shadowgraph of the aerodynamic flow field of an Apollo-Saturn model in a 
wind tunnel, showing turbulence, shock waves, and separated flow. 


Less information is available on fluctuating pressure fields in the hypersonic 
regime (important though they are during atmospheric entry) , other than that 
they tend to cause less structural vibration (refs. 21 and 22). The exception 
is atmospheric entry of low-drag vehicles, where high vibration is observed at 
the lower altitudes (ref. 23). 

Structural protuberances into the aerodynamic flow should also be considered, 
because they can cause added turbulence and fluctuating pressures. These 
protuberances may include ullage, control and retro-rockets, stabilization fins 
electrical cabling and propellant lines or their shrouds, and communication 
antennas . 




4. ANALYTICAL DETERMINATION OF VIBRATION 

Various methods and combination of methods are available for determining 
structural integrity and predicting equipment -vibration requirements, 
especially during the early phases of vehicle development. These methods may 
be divided into the following three categories: (1) classical analysis, 

(2) statistical energy analysis, and (3) extrapolation. 

4.1 CLASSICAL ANALYSIS 

Deterministic or classical analyses are those in which (a) the vehicle structure 
or one of its sections is represented by a mathematical model; (b) the applied 
loading is described by its time history or loading spectrum, and in the case 
of aeroacoustic loading, by its spatial distribution over the exposed surface; 
and (c) the resulting vibration at various locations of the structure is 
calculated by the solution of equations of motion derived from the model. De- 
scriptions of analytical methods may be found in references 24 through 29 . 

4.1.1 Selection of Mathematical Model 

The selection of the mathematical model is influenced by the desired accuracy 
and frequency range, the details to be employed in describing the loading and 
the structure, and the cost of computation and model formulation. The mathe- 
matical model is an important element in vibration calculations. Careful con- 
sideration is usually given to stiffness and mass distributions, and boundary 
considerations in the synthesis of the model. 

Basically, two types of models are available: (a) continuous - or distributed- 

parameter representations and (b) lumped- or discrete-parameter representations. 
Distributed representations are mostly used for fairly simple structures for 
which classical solutions are known, such as domes, conical and cylindrical 
shells, plates, and beams, all with simple boundary conditions (refs. 28 through 
30). In addition, parametric studies of structural vibration often use a 
distributed representation (ref. 31). It is often found, however, that 
composite construction, complex boundaries (including joints and cutouts) , and 
attached masses make a distributed analysis difficult to formulate, often 
inaccurate, and expensive or impossible to solve. Thus a lumped-parameter 
representation is usually preferred to model the structure. In addition, one 
form of lumped-parameter modeling, called finite-element modeling, is gaining 
wide usage, mainly because of its adaptability to matrix notation and high-speed 
digital computation (refs. 26, 32, and 33). Figure 4 shows a typical lumped- 
parameter model of a space-vehicle section. Rigid or lumped masses, whose 
weights are selected on the basis of the mass distribution of the structure, 
are not shown, but are placed at all connecting points of the model. 

The accuracy of the resonant frequency, mode shape, and vibration stress and 
motion calculations will generally be highly dependent upon the number and 
location of these masses and the frequency considered. For the lower-order 
modes, some of the structural details may be compromised, while others may not. 
For example, a structure represented by a lumped-parameter model requires fewer 
lumped masses for the lower-order modes without significantly influencing the 
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Figure 4.-A lumped-parameter model of a section of a space vehicle. 

proper estimation of the resonant frequencies and mode shapes. At locations of 
high vibration stress (and thus, possible failure), however, the local details 
are vital for proper estimation. On the other hand, the structure requires 
more lumped masses for the higher-order modes to preserve resonant frequency 
and mode- shape accuracy, but there is often less likelihood of high stress 
points and thus fewer of these local details to consider. If many masses are 
used for the higher-order modes, the cost of computation is almost always 
unusually large and often overwhelming. Thus, the details of the lumped-para- 
meter model should be selected with great care. 

All real structures vibrate nonlinear ly, so that the transfer function (i.e., 
the ratio of the response to the applied loading) varies with the magnitude of 
the loading. Unfortunately, there have been almost no nonlinear vibration 
analyses of mathematical models that approach the complexity of nearly all space- 
vehicle structures (ref. 34). However, as pointed out by Lyon (ref. 35), non- 
linearities need not be considered unless their effects rival the uncertainty 
of the linear estimate. When sufficiently severe nonlinear behavior is 
expected, it is common practice to assume a linear model and perform the 
classical analysis using the linear values of structural stiffness and damping 
which represent the expected, or the minimum, effect of the nonlinearity on the 
response. The linear response is then scaled, using results of nonlinear 
vibration studies performed on simple models, such as those reviewed in refer- 
ence 34, to provide an approximate or conservative estimate of the nonlinear 
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response. It is assumed hereafter that the nonlinearities are sufficiently 
small to permit the use of linear models. 

4.1.2 Formulation and Solution of Equations of Motion 


Regardless of the specific method of analysis used, such as those described in 
references 24 and 25, for a linear structure the mathematical formulation of 
the problem will lead to a set of equations of the form 

[M]{q> + [C]{q> + [K]{q> = {Q} (1) 

where [M] , [C] , and [K] are the square matrices of mass (or inertia), damping, 
and stiffness coefficients, respectively (and are called the mass, damping, and 
stiffness matrices); and {q} and {Q} are the column matrices of the coordinate 
displacements and the applied forces, respectively. If the damping matrix can 
be diagonalized by the same transformation that uncouples the undamped system, 
classical normal modes exist [i.e., in each normal mode, the various locations 
of the structure vibrate in phase or 180 degrees out of phase (ref. 36)]. A 
well-known special case of classical normal modes is Rayleigh’s proportional 
damping, where the damping matrix [C] is a linear combination of the stiffness 
matrix [K] and/or the mass matrix [M] . Otherwise, nonclassical normal modes 
exist, requiring special treatment (refs. 25 and 37). For classical normal 
modes, the solution of equation (1) will yield the resonant frequencies, mode 
shapes, and the vibration displacement to a given loading, which can then be 
used to determine structural adequacy and the motions and forces applied to 
equipment. For nonclassical normal modes, it is usually assumed that they are 
classical in order to avoid certain mathematical complexities in obtaining a 
nonclassical solution. 

4. 1.2.1 Natural Frequencies and Mode Shapes 

The solution of the homogeneous form of equation (1) without damping (i.e., [C] 

= {Q} = 0) provides the resonant or natural frequencies and mode shapes which 
characterize the mathematical model and thus the structure. There are many 
methods available for calculating resonant frequencies and mode shapes. It is 
common to group them into three categories, depending upon their mathematical 
formulation: (1) energy methods, (2) differential equation methods, and (3) 

integral equation methods. Energy methods include Rayleigh’s method (for the 
fundamental mode) and the Rayleigh-Ritz method (for higher- order modes) . 

Differential equation methods include the Holzer, Myklestad, and Thomson methods. 
Integral equation methods include the collocation, Galerkin and Stodola methods. 
These methods have their relative advantages and limitations as to their accuracy 
and their ease or cost of implementation, which are discussed in references 24, 
25, and 33. 

The principal use of modal data is the determination of the vibration response. 
The resonant frequencies and mode shapes are, however, also useful in selecting 
mounting locations for equipment and guidance sensors, such as gyros and acceler- 
ometers. This is especially important when known critical frequencies will 
cause detrimental effects to the equipment, structure, or control system. 
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4. 1.2. 2 Vibration Response 


The response of a structure at various frequencies, or at desired intervals 
throughout the frequency range of interest, is usually calculated by summing 
the response in each of the orthogonal modes which characterize the structural 
vibration. Since the applied loading is generally random in nature, the 
response is also random. The vibration at location x of the structure is 
therefore expressed in terms of its spectral density function G(x,f). Extending 
the work of Wang and Uhlenbeck (ref. 38), references 25, 34, and 39 show that 
the displacement spectral density for location x at each frequency f due to a 
spatially-distributed applied loading is 
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where the cross -joint acceptance function is given by 
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the frequency response function by 


H ± (f) = [1 - (f/f ± ) 2 + iZ^f/fp" 1 


(2b) 


and the modal or generalized mass by 
M. = J ~ [p (x)h(x) ]dx 

1 r* 1 


(2c) 


and where f . , £., and <J>. (x) are resonant frequency, viscous damping ratio, 

-t th 

and mode shape (at location x) in the i mode and Z direction, respectively; 
H*(f) , the complex conjugate of H^(f); A, the area of the surface exposed to 

the fluctuating pressure field; [p(x) h(x)], the surface density (i.e., mass 
per unit area) of the structure at location x; S, total area of structure; 
and G (f) , the reference fluetuating-pressure spectrum (often the spatial 

average over the exposed surface). The cross-spectral density function of 
the fluctuating pressures for locations £ and £ 1 , of the exposed surface is 


G (C,5',f) = C (?,5\f) - iQ (5,c,f) 

lr r r 

where the cospectrum and the quadspectrum are 


(3) 


c (5,?',f) = 2 f [R (C,5',x) +R (5,C',-t)] cos 2irfxdx 

P •/ P P 


(3a) 


20 



Q p (5,5',f) 


(3b) 


= 2 
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[R p a,c’,T) - 


R^(£,£ t ,-t)] sin 27rfxdT 


respectively (ref. 40). 


Thus, the pressure cross- correlation R^ (£,£ f ,x) must 


be measured or assumed as a function of the time delay t for every combination 
of location points £ and £ ! . References 7, 34, 39, 41, and 42 describe and 
discuss several distributions of cross -correlation functions that might be 
assumed for various aeroacoustic-noise loadings. Figure 5 shows the measured 
cross-correlation coefficients, R^ (£ * ,t) / a p (£) a p (? T ) , in various frequency 

bands for a pair of points during flight, where a^(0 is the rms pressure at 

location £. 


Substitution of equation (3) into equation (2) will determine the displace- 
ment spectral density of the vibration response. Computer programs have been 
developed recently to calculate this response (refs. 43 through 46). The 
acceleration spectral density, usually preferred in specifying design and test 

4 

requirements for equipment, can be found from G..(x,f) ~ ( 2 tt f ) G (x,f). Spectral 

w w 

density equations for moments and forces, from which the stresses throughout 
the structure can be calculated, are generally more complicated than equation 
(2) for displacement and appear in the appendix. 



Considerable simplification of equation (2) can be achieved if the viscous 
damping ratio in each mode is small, and the resonant frequencies f^ are 

well separated. Then at the resonant frequencies 


G w (x, f i) ~ A 2 G pr (f ± )(j)^(x)Q^j^(f i )/(2irf i ) 4 M^ 


(4) 


1 th 

where = resonant magnification or quality factor in the i mode. 


The analyst is confronted with two major problems in achieving an accurate 
prediction using equation (2) or (4) . The first problem is the proper 
determination of the pressure cross -spectral density G (£,£',f) or cross- 


correlation R^(£;,£ t ,t), used in calculating the joint acceptance. If the 

cross -correlation is assumed, the assumption may be invalid. If it is to 
be measured, it will be necessary to use a large number of wideband measurement 
channels during wind tunnel or flight test. Also, there will be an additional 
cost to provide accurate relative phase characteristics throughout the data 
acquisition and reduction process. Wind tunnel measurements may require the 
use of a low-noise wind tunnel. Otherwise, high tunnel noise may mask more 
moderate aerodynamic noise over the wind tunnel model. Also, the microphone 
size must be scaled. Since flight testing usually occurs late in the vehicle- 
development program and to avoid the cost of telemetered flight data, the 
cross-correlation is usually assumed, often inaccurately. 


The second problem is the proper determination of the viscous damping ratios 

in the various modes. In equations (2) and (4), viscous damping is assumed 

because of its compatibility with the assumption of linear vibration and the 

subsequent ease in solving equation (1). Unfortunately, the damping of real 

aerospace structures is not viscous, but usually occurs in one or more of 

the following three forms: (1) material damping, (2) friction damping, and 

(3) acoustic radiation. Material damping is the inelastic behavior of the 

structure caused by internal friction during microscopic slip between interfaces 

within the nonhomogeneous material. Lazan has shown experimentally that the 

damping in a mode of vibration is dependent on the stress distribution through 

the material, i.e., on the maximum stress, as shown in figure 6, and the mode 

shape (f) . (x) (refs. 47 and 48). If the maximum vibration stress s is less 
i max 

than the endurance limit stress s of the material (s <0.7 s in fig. 6), 

e max e 

the damping changes only slightly with stress. If the maximum stress (from 
the contribution of all modes) exceeds the endurance limit, the damping 
increases significantly for most materials; but how it is distributed in the 
various modes is presently unknown. However, because of the increased damping, 
the life of the structure is prolonged. Friction damping is energy dissipa- 
tion caused by slip or sliding between mated surfaces. Slip is relative motion 
within a region of the mated surface, whereas sliding occurs over the entire 
surface. Since most aerospace structures are comprised of parts that are 
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Figure B-Material-damping properties for various aerospace structures. 


bolted or riveted together, friction damping is usually the dominant source 
of damping. The friction force per unit area during slip or sliding is 
approximately constant and is controlled by the pressure normal to the mated 
surfaces (refs. 47, 49, and 50). Because the friction force is not 
proportional to velocity, as is the viscous-damping force, a significant 
nonlinearity exists. Since it is approximately constant, friction damping 
is unfortunately less effective for controlling high vibration. Also, the 
static friction must be overcome for slip or sliding to occur. Thus, there 
is little or no friction damping at low vibration. The proper estimation 
of friction damping and its distribution in the various modes is very difficult. 
This is because of the distribution of the normal pressure over the mated 
surfaces between the fasteners, the variations between the forces exerted by 
the fasteners (due to the lack of control during manufacturing and fabrication) , 
and the spatial variations in the static and kinetic coefficients of friction. 

Acoustic radiation (i.e. , air damping) is the generation of acoustic waves by 
the vibration of the structure and its propagation to other structures and the 
surrounding space. Radiation can be an appreciable form of damping for 
structures with low surface densities, such as panels, unless the structure 
vibrates at high altitude (where the radiation is reduced) or the surrounding 
space is reverberant. For simple cases, such as acoustic radiation from a 
rectangular panel into a free field, or a reverberant space with known surface 
characteristics, radiation damping in the various modes can be properly 
estimated with relative ease (refs. 51 through 56). In the more complex cases, 
the estimation can be quite involved. 



There are also other forms of damping. One of them is viscoelastic, which 
damps by the inelastic behavior of certain rubber- like materials caused by 
nonlinear distortion of long- chain molecules. These materials are sometimes 
bonded as a layer or sheet to the surface of structural panels (often in 
combination with foil backing as a tape), or are sandwiched between adjoining 
structures in the joints or interfaces. References 47 and 57 may be used to 
estimate the contribution of viscoelastic damping. 

Adjacent structures that have small spaces between facing surfaces, as might 
be found in joints that connect panels with stiffeners or other panels, can 
exhibit damping in the form of air “pumping. 5 5 (This dissipation is in 
addition to the friction, radiation, and viscoelastic damping discussed 
previously#) The contribution of air-pumping is described by Maidanik (ref. 

58) . 

From this discussion it should be obvious that the proper estimation of the 
damping in the various modes will be difficult or impossible, except for the 
simplest structural configurations. Thus, to determine damping, there is 
usually great reliance on data from previous tests of similar structures. 

Even then, the data must be extrapolated to account for the differences between 
the structural configurations and the test and flight conditions. A conserva- 
tive value is usually selected for the “equivalent viscous 5 * damping ratio 
to be used in equation (2) or (4) . Sometimes an iteration procedure is used 
to modify the damping selection after the vibration displacement and stress 
are calculated in the various modes. Nonlinear representation of damping in 
complex structures is considered to be beyond the present state of the art. 

Vibration response is often caused by excitation from more than one source 
(e.g., aeroacoustic noise applied to the external vehicle surface, and 
mechanically transmitted vibration direct from the engines) . During a parti- 
cular phase of the mission, if only one source is dominant, as illustrated in 
figure 1, then only the response to that source need be considered. On the 
other hand, if two or more sources are effective during the period, multiple- 
input analysis may be used, as illustrated in figure 7. Reference 40 shows 
that the response spectral density for location x at each frequency f due to 
multiple loading is 

G (x,f) ^ ^ G IK^ • ^ 

1=1 K=1 


where H (f) is the transfer function between the response (displacement, 

JLX 

acceleration, moment, stress, etc.) at x and the applied loading (force, 
pressure, acceleration, etc.) at I; and G_^(f ) , the cross -spectral density 

between sources I and K. The units of H^(f) will, of course, differ for 

each type of loading and response. 
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Figure 7-Diagrammatic for determining vibration response to multiple sources. 

Usually, the sources are independent of each other. Then equation (5) reduces 
to 

A ? 

G (x,f) = V IS Ix < f ) I GjCf) (5a) 

1=1 

Values of the transfer function may be determined analytically or experimen- 
tally. As shovm in references 34, 39, and 59, the rms response a (really the 
standard deviation) can be obtained from the variance 



(5b) 


When the applied loading is sinusoidal, which it is occasionally, the displace 1 
ment at location x and frequency f of the sinusoid is 


w(x,t) = w (x) cos (27rft + 0) = H (f) F (5) cos 2irft (6) 

O “lx o 

where w q and 0 are the amplitude and phase of the response displacement, 
respectively; and F q (£)> the amplitude of the applied loading at location 
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4. 1*2.3 Failure Modes 


Despite the fact that most space vehicles undergo relatively short exposure 
time to high vibration during flight, the most common cause of structural 
vibration failures is fatigue. The major factor is the relatively high fre- 
quency of the vibration. A great deal of experimental data are available on 
fatigue under conventional single- stress amplitude single- frequency sinusoidal 
stresses, and combined static and sinusoidal stresses (refs. 60 through 63). 
Some fatigue data are available for multistress amplitude single-frequency 
( € ‘programmed* 5 ) stresses, but little for multifrequency or random stresses 
(ref. 64). In most cases, therefore, sinusoidal fatigue data must be 
extrapolated for use in estimating damage and possible failure under random 
vibration stress, as found in space vehicles. 

Sinusoidal fatigue data are conventionally plotted as s^ vs N_^ curves , as 
shown in figure 8, where s^, s^, and s^ are stress amplitude, ultimate ten- 
sile stress, and endurance limit, respectively; and n_^ and ISL, the number of 
applied cycles and number of cycles to failure at stress amplitude s^, respec- 
tively. The s^ vs I\L curve is dependent on the static stress, temperature, 

and frequency, since they affect the fatigue properties. If the stress 
contribution from one of the modes is much greater than the contribution from 
all other modes, the resulting narrow-band random stresses are quite similar 
to multistress amplitude single-frequency stresses. In these cases, it is 
possible to make two assumptions which are useful in calculating fatigue 



NUMBER OF APPLIED CYCLES (n ; ) 
NUMBER OF CYCLES TO FAILURE (Nj) 


Figure 8 -Fatigue curve for a typical aerospace material under sinusoidal loading. 
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damage. First, it can be assumed that the damage D accumulates as in the 
Palmgren-Langer-Miner hypothesis, which states that D = ^ C^/N^) • Using 


this hypothesis, fatigue failure will occur when D = 1. Second, it can 
.be assumed that the stress peaks vary in accordance with the Rayleigh 
distribution: 


p p (s ± ) = exp [-s 2 /2 cr 2 ] (7) 

where a is the rms stress. All of the various damage -accumulation hypotheses 
s 

produce errors, the magnitude of which varies with the material and the 
sequence of the stress amplitudes applied. The Palmgren-Langer -Miner hypothesis 
is commonly used because it is no worse than the others, such as the Corten- 
Dolan hypothesis (ref. 65), and it is easy to apply. The actual distribution 
of random peaks generally approximates the Rayleigh up to fairly high stresses 
(s^ < 3a g ) , beyond which the Rayleigh assumption is almost always conservative. 

Applying a modification to the analysis first performed by Miles (ref. 66), 
the damage under narrow-band random stresses is 

D = [f n T/( Sl /a s ) 2 ^][2^ r(q+1)][I(u u ,q) - I(u e ,q)] (8) 

where f^ is the resonant frequency of the mode which has the dominant stress 

contribution; T, the exposure time of the applied loading; q = 3C a/ 2 (defined 
in the next paragraph); T(q+1), the complete gamma function (ref. 67); and I(u,q) , 
the ratio of the incomplete to the complete gamma function shown in figures 9 

and 10 (ref. 68), u = (q+1 ) ^ ^ s^/2 , and u = (q+1) ^ ^ s^/2a^ 

u us* e es 

Crandall (ref. 59) has shown that the damage variability due to this 

randomness is usually small. 

Material fatigue properties s^ and a are shown in figure 8 as the ordinate 

intercept and slope of the s ± vs N curve. The sine- to-random fatigue 

conversion factor necessary to provide random fatigue failure when D=1 , 

varies with the material. Its value has been observed to vary over a wide 
range (1/2) <J^<2, and is the subject of some disagreement (refs. 66 and 69). 

If more than one mode contributes to the stress, the distribution of peaks is 
a combination of the Rayleigh and the Gaussian distributions (ref. 70). Unfor- 
tunately, no widely acceptable damage theory has been developed for this case 
(ref. 64). However, the use of equation (8) will provide a conservative 
estimate of damage (ref. 71). 

Another possible failure mechanism is the exceedance of the ultimate dynamic 
stress of the material by a single stress peak. This is a special case of 
the random-vibration problem, commonly called the threshold-crossing or first- 
passage problem, and specifically the probability P t ( 3>3 q ) that a single 
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stress peak s. exceeds or equals the threshold (or ultimate dynamic stress) 

* 1 * 
s during the exposure time T, where 3 = s./o , 3 = s /a , and a is the rms 

u is ous s 

stress. If the stress contribution in one of the modes is much greater than 
the contribution from the other modes, the probability may be approximated by 

P T (S ^ 3 C ) = 1 - exp [-27ra o f n Q n T] (9) 

where f^ and are the resonant frequency and resonant magnification of the 

dominant mode, respectively, and a is as given in figure 11 (refs. 72 and 

73) , in which 3 = 3 + should be used when failure may occur in tension only 
(i.e. , a one-sided threshold), and 3 = |^| should be used when failure may 
occur in either tension or compression (i.e., a two-sided threshold). It is 

* 

important to note that the ultimate dynamic stress usually is greater than 

the ultimate static stress s^, since most materials can withstand a higher 

stress when the duration of application is short (ref. 74). Unfortunately, if 
more than one mode contributes substantially to the stress, no experimental 
data are presently available on the probability of threshold exceedance in a 
given time period. 

There are other possible failure mechanisms under vibration, including colli- 
sion between adjacent structures, wear, and operational drift. Collision 
between adjacent structures is another case of the threshold-crossing problem 
and is usually only of concern to elements within equipment rather than to 
primary or secondary vehicle structure. Wear is seldom a failure mechanism in 
space vehicle structures and operational drift is usually associated with 
guidance equipment (refs. 75 and 76). Thus, neither will be discussed in this 
report. 

To determine the adequacy of the structure for all potential failure modes, 

the effects of vibration must be assessed in combination with other loads 

existing simultaneously and sequentially. As an example, high vibration 

may be expected near the q period of flight. During this period, high 

max 

static stresses may be induced by the sustained acceleration and differential 
pressure resulting from venting lag. Since all of these act simultaneously, 
their combined effect must be considered. 

In certain cases, the combined effects may be readily determined. For 
example, s_^ vs N. curves are available for many common metals under combined 

static and sinusoidal stresses and various temperatures, which may be used in 

figure 8 and equation (8) when one mode dominates the random response. Also, 

the probability of exceeding a threshold under combined static and random 

stresses, which may be calculated when one mode dominates by substituting 

3 = (s. + s ) /a into figure 11, where s is the static stress and the 
i — st s st 

+ sign indicates that the static acceleration causes a tensile or compressive 
stress, respectively. If the static stress is large enough, the threshold 
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exceedance is predominantly one-sided, so that $ = should be used. Unfor- 
tunately, in many cases, state-of-the-art advances are required before com- 
bined effects can be predicted (ref. 77). 

In the preceding discussion, it is assumed that single-point structural 
failures must be avoided. However, under certain circumstances, structural 
redundancy allows one or more of these failures without loss of the complete 
structure. These circumstances are described in reference 78. 

4.1.3 Assessment 


Classical analysis can be the most accurate method of predicting vibration 
stresses and motions in the lower- frequency range. Typically a rigorous 
analysis is restricted to that range which encompasses the lower 50 modes of 
the structure or section. The lack of structural detail in the model of 
simplifications made in describing the loading often further reduces the use- 
ful frequency range. 

4.2 STATISTICAL- ENERGY ANALYSIS 

Acoustic and aerodynamic noise cause the vehicle structure to vibrate in many 
modes over a broad frequency range. As noted in Section 4.1, classical 
methods often provide results which are inaccurate and/or expensive in the 
higher-order modes. An alternate approach, statistical energy analysis (SEA) , 
has been developed by Lyon and his associates (refs. 79 through 84) to esti- 
mate the vibration of complex structure subjected to random loading at high 
frequencies. SEA can be used at the high frequencies to provide upper-bound 
and broad-brush estimates of vibration response and transmission with few 
calculations (relative to those required for classical methods) , using only 
gross structural properties. However, SEA and classical methods both suffer 
from the lack of accurate information on the structural damping (material and 
friction) of the various modes. 

4.2.1 Gener al Sta tist i cal-En ergy Analysis Formulation 

The primary variables in statistical energy analysis (SEA) are power flow and 
modal energy (i.e. , the vibration energy per mode). A fundamental principle 
of SEA states that under wideband random loading: 

• The time -average power flow between two coupled dynamic 
systems is proportional to the differences in the time- 
average modal energies of the systems. 

• The power flows from the high-modal- energy system to 
the low-energy system (ref. 52). 
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When system 1 is excited by wideband random loading and system 2 is excited 
only through coupling, the relationship between the time-average modal ener- 
gies 0^ and of the two systems is 

9 2 /® i = n 21^^ n 21 + n 2^ 0°) 


where r\^ an h are t ^ ie damping of systems 2 and the coupling loss factor 

from system 2 to system 1, respectively. Since > 0 and 1121 > 0» it 
results that 0^ < 0^* From this result, a rule can be formulated for esti- 
mating an upper bound on the vibration of unexcited structures that are 
coupled in cascade to the directly excited structure; namely, the time- 
average modal energy of these unexcited structures cannot exceed the time- 
average modal energy of the directly excited structure. 


SEA is used to predict the space- average response of a system, rather than 
the response at particular locations. The modal energies are simply related 
to common dynamic parameters. For a flat plate, the modal energy 0^ is 

related to the space-average acceleration spectral density <G^.(f)> by 


0 1 = [ph 2 c^/3 2 (27rf) 2 ] < G^ (f) > 


(10a) 


where p, h, and c' are mass density, thickness, and speed of longitudinal 

Li 

waves of the plate, respectively. For a reverberant acoustic field, the modal 
energy 0^ is related to the space-average pressure spectral density <G^(f)> 

by 


0 = [c /p (2irf) ] < G (f ) > 

ac a a p 


(10b) 


where p^ and c^ are mass density and speed of sound in the acoustic medium 
(i.e., air), respectively (ref. 81). 

4.2.2 Applicabili ty to Simple Vehicle Structures 

SEA calculations are quite simple for many vehicle configurations. Figure 12 
shows an idealized representation of a simple structural configuration often 
found in space vehicles. Acoustic noise from the rocket engine causes the 
external structural section to vibrate. Because of the spatial distribution 
of the noise sources and the sound- reflection effects of the launch pad or 
test stand and the surrounding terrain (discussed in Sec. 3.1), the acoustic 
field adjacent to the external structure may be represented by a reverberant 
field with a space-average pressure spectral density <G (f)>^. The external 

, represented here as a flat plate, responds as a reverberant 


structure 




INTERNAL 

STRUCTURE 


Figure 12 -Idealized representation of a simple space-vehicle structure used for predicting the 
space-average response of external and internal structure at higher frequencies by 
statistical-energy analysis. 


vibration field with a space-average acceleration spectral density <G^.(f)>^ 

and transmits vibration to the internal structure (or equipment) * which has 

a space-average acceleration spectral density <G..( f)>.. (Since most external 

w 1 

vehicle sections are cylindrical shells rather than flat plates, curvature 

effects must usually be considered below the ring frequency f = c'/27ra, 

r Li 

where a is the radius of the cylinder,) 


If the vibration power transmitted from external to internal structure is 
small, compared with the total power dissipated and radiated by the external 
structure, the last three equations may be used to determine the resonant- 
motion acoustic acceptance (or vibroacoustic transfer function) of the 
external structure: 


< V f > 

< G p <£> > 5 


<p e h e>' 
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02 ) 


The upper bound of the acceptance is obtained by assuming r| >>n (i.e., the 

ea e 

acoustic radiation damping exceeds the structural damping of the external 
structure). This upper bound is shown in figure 13 for an 0.1 25 -in. thick 
aluminum plate. The coupling loss factor can be calculated from 
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Figure 13. - Space-average acoustic acceptance for an idealized representation 
of space-vehicle external structure. 


n = R j/Zirfm, where R , is the radiation resistance, and m, the mass of 
ea rad 9 rad * * 

the plate (ref. 82). Figure 13 also shows the resonant acceptance of an 

-2 

assumed structural damping of ri e - 10 and the calculated coupling loss 

factor n . In this case, the acceptance peaks at the acoustic-structure 

ea 1/22 

coincidence frequency f =3 c /irh c_! . 

^ J c a e L 

The of f -resonance behavior of the external structure is controlled by its 
mass. The mass- controlled acoustic acceptance for the flat plate is 
2 

2 /(P e ^e^ an( ^ iS a ^ so s ^ own i- n figure 13. The total acceptance is the sum of 
the resonant and mass- controlled acceptances. 


The vibration- transfer function from external to internal structure, obtained 
from equations (10) and (10a), is 


* V £ > V f > % ■ <»l h 2'»e h e> 


n ie /(n l e + V 


( 12 ) 


when the internal structure is represented by a box comprised of six plates, 
as shown in figure 12, and where p^, tu, and are mass density, thickness, 

and structural damping of these plates, respectively. If the internal 
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structure is attached to the external structure by four rigid studs, the 


coupling loss factor is n ie - ztl ^ c L 


, 1 / 2 , 


= 2h^c£ /3 '“fA^, where A^ and c^ are total sur- 
i x 

face area and speed of longitudinal waves of the internal structure, respec- 
tively (ref. 83). The vibration transfer function is shown in figure 14 


for p i = p e’ 


h i - h e’ 


= 10 


-2 


and the calculated value of using 


A. = 12 ft 
x 


This example demonstrates the application of SEA relationships 
and calculations to simple vehicle conf igurations . 


4.2.3 Applicability to Complex Vehicle Structures 

SEA techniques have also been applied to large complex vehicle configurations. 
Sections of the Saturn V launch vehicle have been analyzed to predict vibra- 
tion for a variety of fluctuating pressure fields: acoustic noise at liftoff, 

subsonic and supersonic boundary -layer turbulence, shock-induced separation, 
and disturbed flow at various times during the ascent phase of the mission 
(ref. 84). The launch vehicle was represented as an orthotropic cylindrical 
shell with ring frames and longitudinal stiffeners. The vibration acceleration 
and strain were calculated for each loading condition. It was observed that 
the time- average power from the aerodynamic-pressure field was relatively 
insensitive to the details of the correlation field used to represent the 
pressure source. 



Figure 14.- Vibration-transfer function between external and internal structure for an idealized 
representation of a space vehicle. 
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Acoustical and mechanical vibration transmission has been investigated for a 
simplified physical model of the 0G0 spacecraft (ref. 81). The external 
acoustic field was again represented as a reverberant field, which excited 
vibration of the payload shroud. The vibration was transmitted to the space- 
craft by way of two paths. The acoustic path was comprised of the shroud, 
which radiates acoustic noise into the interior space; the internal reverber- 
ant acoustic field; and the spacecraft. The mechanical path was comprised of 
the shroud, a ring frame, four mounting trusses, and the spacecraft. The 
analysis showed that the acoustic path was more efficient at low frequencies, 
and the mechanical path at high frequencies. An experimental study con- 
firmed this conclusion. 

SEA calculations for both the Saturn V and 0G0 cases were longer than those 
shown previously for simple structures. However, they still represented a 
small fraction of the calculations that would have been required if classical 
methods had been used to provide similar estimates. 

4.2.4 Assessment 


Although still in its early development, SEA is already a valuable tool for 
providing physical interpretation and rough estimates of vehicle vibration 
for frequencies sufficiently above the lower resonant frequencies. Future 
developments are likely to produce improved accuracy and/or greater simplifi- 
cations for complex structures. However, it is unlikely that the simplicity 
of the extrapolation methods discussed in Section 4.3 will be challenged. 

4.3 EXTRAPOLATION METHODS 

Extrapolation or empirical methods utilize experimental data obtained on a 
vehicle, usually called a reference vehicle, to predict the vibration of a 
new vehicle. Certain assumptions must be made for the extrapolation method 
to be used. Often these assumptions are only partially valid, so that errors 
may be expected. Also, presently available extrapolation methods are gener- 
ally used to predict vibration motion (usually in terms of acceleration vs 
frequency for sinusoidal excitation, or acceleration spectral density vs 
frequency for random excitation) rather than vibration stress. Thus, these 
methods are used only to determine vibration requirements for equipment, 
rather than to determine the integrity of structure. The methods may be 
further divided into two subcategories: frequency response methods and scal- 
ing methods. The relative merits and limitations of these methods, as well 
as the classical and SEA analyses previously described, are discussed in 
references 85 and 86 in detail. 

4.3.1 Frequency- Response Methods 


Frequency response methods use measurements obtained on the reference vehicle 
to determine C^(f) , the ratio of the resulting vibration to the magnitude of 

the fundamental source, as a function of frequency or bandwidth. For all these 
methods, the fundamental source is acoustic noise or aerodynamic noise applied 
to the vehicle surface, or both. 
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All frequency response methods use statistical analysis to determine a suf- 
ficiently conservative value C^Cf) to be used as vibration criteria for 
the new vehicle: 


G.. (f) = C (f) G (f) (13) 

wn r crpn 

where G„^(f) is the acceleration spectral density of the structural vibration, 

G (f) the pressure spectral density of the aeroacoustic noise as measured on 

the external structure, and the subscript cr denoting criteria based on the 
reference vehicle. 

4. 3. 1.1 Mahaf fey -Smith Method 

Mahaffey and Smith (M-S) (ref. 87) were the first to apply this method, using 
data obtained on B-52 and B-58 aircraft. First, acoustic measurements were 
made at a sufficient number of locations to define adequately the distribution 
of fluctuating pressures over the vehicle surface, and were analyzed in octave 
bands to determine their frequency content. Figure 15 shows the resulting 
distribution of sound-pressure levels for the octave band from 300 to 600 Hz 
over the surface of the B-58 during ground runup of its turbojet engines. 
Vibration measurements were then made at various locations throughout the 
vehicle, including structure adjacent to points of equipment attachment, and 
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Figure 15rSound- pressure-level contours for the B-58 in the 300- to 600-Hz octave band 
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were similarly analyzed in octave bands. These data were next used to plot 
vibration magnitude vs sound pressure level, as measured on the vehicle sur- 
face closest to the vibration measurement, in each of the octave bands 
throughout the frequency range of interest, (i.e., a lower limit of 37.5 Hz 
and an upper limit of 2.4 kHz). Figure 16 shows this plot for the 300- to 
600-Hz octave band. 

As seen from the first ordinate in figure 16, the vibration data are plotted 
in terms of peak rather than rms values, even though acoustically-induced 
vibration from jet engines is invariably random. M-S were required to present 
their data in peak values to derive sinusoidal design and test specifications, 
as required on the B-58 by the Air Force. However, M-S stated that generally 
the ratio of the peak value to the rms value was about 3.3. Because the 
vibration varied considerably throughout the vehicle, even for the same 
sound pressure level, the remainder of the octave band plot showed a con- 
siderable dispersion of data (fig. 16). This would be expected since there 
was no effort made to categorize the vibration data according to sections of 
the vehicle (called zones), the type of structure, or direction monitored. 
Thus, higher vibration would be expected for external panels, wing locations, 
unloaded structure, and directions normal to the nearest vehicle surface, 
than for internal structure, fuselage locations, equipment- loaded structure, 
and directions tangent to the nearest vehicle surface. M-S handled this 
dispersion statistically, using the plotted data to establish the regression 
line, and, assuming a log-normal distribution, the standard deviation. 



Figure 1 6.— Vibration-prediction curves for Mahaffey-Smith and Brust-Himelblau methods in the 
300- to 600-Hz octave band 
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Figure 16 shows the regression and various confidence lines for the 300- to 
600-Hz octave band. It is also noted that for any of these lines, the vibra- 
tion increases by a factor of about 4 (rather than 10) for a 20-dB increase 
in sound-pressure level. This demonstrates that, statistically speaking, the 
B-58 structure is nonlinear. 

To use the M-S method for a new vehicle, it is necessary to assume, among 
other things, that (a) the vibration is primarily acoustically induced, 

(b) the new vehicle and the B-58 are dynamically similar statistically, (c) 
the lack of data categorization has no particular disadvantage, and (d) the 
desired unit of vibration magnitude is the peak acceleration in an octave 
band. Then the acoustic noise spectra at the locations of interest on the 
surface of the new vehicle must be established, using the references noted 
in Section 3. 

4. 3. 1.2 Brust-Himelblau Method 

Brust and Himelblau (B-H) (ref. 88) reconverted the M-S vibration measure- 
ments from peak values to the original rms values (a factor of 3.3), then 
squared them and divided by the octave bandwidths to express the random 
vibration in more suitable units, the (mean square) acceleration spectral 
density (ASD) averaged over the octave band. The second ordinate of figure 16 
shows this conversion. Following a M-S suggestion, B-H used the upper 60% 
confidence line to derive vibration design and test requirements for equip- 
ment locations. However, it is common knowledge that random vibration 
measured in narrow bands varies considerably from that averaged over an 
octave band. 

Using the narrow-band, one- third-octave band, and octave-bandd analyses of 
random vibration measured on the Titan I missile, shown in figure 17, B-H 
concluded that the ratio of narrow-band ASD to octave-band ASD was nearly 
always less than a factor of 5. Thus, if vibration requirements are to be 
derived by enveloping the narrow-band peaks, the octave-band-averaged ASD 
from figure 17 must be multiplied by five to provide the necessary 
correction. 

To use the B-H method, it is necessary to make the assumptions used in the 
M-S method, except that the desired unit of vibration magnitude is either 
the octave -averaged or the narrow-band ASD. For the narrow-band ASD, the 
assumption relating the narrow-band ASD to octave-band ASD- must be acceptable. 

4. 3. 1.3 Eldred, Roberts, and White Method No. 1 

Eldred, Roberts, and White (E-R-W) (ref. 4) analyzed vibration and acoustic 
measurements made on the Snark missile and plotted them in a manner almost 
identical to that of M-S. The random vibration was, however, plotted in 
terms of rms acceleration rather than peak acceleration in each octave band. 
Because fewer measurements were used, only the regression line was drawn. 
Converting to compatible units (a factor of 3.3), the E-R-W regression line 
was always found to be higher than the M-S regression line and usually varied 
between the M-S upper 50% and upper 95% confidence lines. Also, above 150 Hz, 
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Figure 17. - Comparison of narrow-band, one-third octave-band, and octave-band spectral 
analyses of the same random-vibration measurement. 

the rms vibration increased by a factor of 10 for a 20-dB increase in sound- 
pressure level, demonstrating that the Snark structure was linear on a 
statistical basis. 

To use the E-R-W method, most of the assumptions mentioned under the M-S 
method must apply. In addition, the E-R-W restricts its use to the regres- 
sion line. The E-R-W method could be revised to more suitable random- 
vibration units (i.e., ASD) , using a conversion similar to that of B-H on 
the M-S method. Also, a narrow-band correction could be applied to the octave- 
band vibration, as mentioned by B-H. 

4. 3. 1.4 Eldred, Roberts, and White Method No. 2 

In addition to the method just described, Eldred, Roberts, and White devised 
a second method (refs. 4 and 89) based on the response of a simple mechanical 
oscillator (i.e., a single-degree-of-freedom system) to a random force applied 
to the mass 

4“* f n Q V f n )/2m2 (14) 

where is the rms -acceleration response; f^, Q, and m are resonant frequency, 
resonant magnification or quality factor, and mass of the oscillator. 


41 


respectively; and G (f ) is the force-spectral density of the random force at a 

r n 

frequency corresponding to the natural frequency of the oscillator. Great 
care must be exercised in the use of equation (14) because the rms accelera- 
tion is infinite when the random force has a white spectrum. However, as 
shown in reference 54, equation (14) generally provides a reasonable approxi- 
mation when the random force has a band -limited white spectrum and the 
oscillator has a Q > 20. 


To apply this equation to a multimodal structure which is typical of space 
vehicles, E-R-W suggested that the right-hand side of equation (14) be 

2 

multiplied by a multimodal-response factor 6 . To determine the value of this 
factor to be used, E-R-W examined vibration and acoustic data obtained on 
eight vehicles: the Corporal, Vanguard, Thor, Titan I, Talos , Polaris, Atlas, 

and Terrier. They assumed that (1) the effective mass m was the gross mass of 
the entire vehicle, (2) the effective random force G^f^) was the product of 

the average pressure spectral density (over the entire vehicle) times the 
square of the entire surface area, and (3) Q = 15. 


2 

Figure 18 shows 6 as a function of the normalized wave number (irfD/c ) where 

a 

D is the vehicle diameter and c , the speed of sound in air (c =1130 ft /sec 

a 2 a 

at sea level and at 70 F) . As seen from figure 18, 6 » 1 . 6 on the average 

and did not exceed a factor of 50 for the vehicles surveyed. To use this 


method, it is necessary to assume that (a) the vibration is primarily 



Figure 18. - Contribution to the rms acceleration from multimodal response as measured on 
several missiles. 
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acoustically induced, (b) the lack of data categorization is not a particular 
disadvantage, and (c) the desired unit-of -vibration magnitude is the rms 
acceleration over the entire frequency range of interest. 

Since the eight reference vehicles have a wide dispersion of dynamic 
properties, the assumption of dynamic similarity between the new vehicle and 
the reference vehicles can generally be considered satisfied. However, on a 

2 2 
theoretical basis, 6 should never be less than unity. Since 6 <1 was cal- 
culated from many data points, it must be assumed that the error was caused 
by the assumptions of the effective mass, random force Q, or other factors 
discussed in reference 4. 

4. 3. 1.5 Curtis Method 

Curtis (ref. 90) observed that the vibration (caused by flow-induced aero- 
dynamic noise) measured on aircraft during high-speed flight was related to 
the aerodynamic pressure q. The values of q studied varied from 90 to 1760 
psf. He performed a statistical analysis of vibration data obtained from 
F8U, B-59, F-101 and F-106 aircraft and concluded that the data approximated 
a Rayleigh distribution. Curtis established that 

G,. (f) = (q/2130) 2 G..(f) (15) 

wn n wr 

where G..(f) is the acceleration spectral density of the structural vibration 
w 

and n and r the subscripts denoting new and reference vehicles, respectively. 


For structure adjacent to internally located equipment, it was determined 

2 

that the mean reference spectral density G ^ r (f) = 0.006 g / Hz, and that 98% 

of the reference narrow-band spectral peaks (based on a Rayleigh distribution) 

2 

were less than the recommended criteria G.. (f) =0.11 g / Hz, 20 Hz<f<2 kHz. 

far y cr e * 

For structure adjacent to equipment exposed directly to the aerodynamic flow 

2 

(such as IR sensors and antennas), G„^(f) =0.011 g / Hz, 20 Hz<f<150 Hz, and 

0.020 g^/Hz, 150 Hz<f<2 kHz. The 98-percentile G ^ r (f) cr “ 0.130 g^Hz, 

20 Hz<f<150 Hz, and 0.235 g^/Hz, 50 Hz<f<2 kHz. To use the Curtis method, it 

is necessary to assume that (1) the vibration is induced primarily by aero- 

dynamic noise, (2) the acceleration spectral density of the random vibration 

2 

is directly proportional to q , and (3) the new vehicle is dynamically similar 

to the four reference aircraft on a statistical basis. Only the aerodynamic 

pressure q^ of the new vehicle need then be established. 

4. 3. 1.6 Franken Method 

Franken (ref. 91) performed a similar study of vibration and acoustic meas- 
urements of Jupiter and Titan I external structures (stiffened cylindrical 
shells) exposed to acoustic noise from rocket engines during static firing. 



analyzing the data in one- third- octave bandwidths. Vibration measurements 
were made in the radial direction only (i.e., normal to the surface). For 
each band, the rms vibration was multiplied by the average weight -per-unit 
area of the structure (called the average surface density) and then divided 
by the rms-acoustic noise, which is one form of a frequency -response 
function (FRF). The use of the weight-per-unit area is related to * ‘mass- 
law 5 * relationships in classical acoustics. 


The vibration of a cylindrical shell can be related to the ratio of the 
excitation frequency to the isotropic ring frequency, that is, (f/f^) = 

(fuD/c^) , where D is the vehicle diameter and c£, the speed of longitudinal 

waves in the material (approximately 17 000 to 18 000 ft/sec for most aero- 
space materials, and 40 000 to 50 000 ft/sec for beryllium and certain boron 
and carbon composites). Taking this as a cue, Franken plotted his FRF vs the 
product fD to obtain the dashed curve shown in figure 19. A maximum scatter 
of 6 dB is observed in this figure. Also, the maximum FRF occurs at fD = 
6000, corresponding to the ring frequency. To use the Franken method for a 
new vehicle, it is necessary to assume that (a) the vibration is primarily 
acoustically induced, (b) the new vehicle and the Jupiter and Titan I 
external structures are dynamically similar within a 6-dB scatter, and (c) 
only the radial direction is of interest. Then the acoustic noise spectra 
at the locations of interest on the surface of the new vehicle, as well as 
its diameter, must be established. However, if vibration requirements are 
to be derived by enveloping the narrow-band peaks, then the one -third -octave 
band-averaged ASD must be multiplied by a factor to provide this correction. 
From the comparisons shown in figure 17, B-H suggested a correction factor 
of 2.5. 



Figure 19. - Vibration-prediction curves for Franken and Winter methods. 
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4. 3. 1.7 Winter Method No. 1 


Winter (ref. 92), who assisted in the development of the Franken method, has 
modified the method after adding Minuteman, Skybolt, and Genie data to the 
previously used Jupiter and Titan I data. As shown in figure 19 as the solid 
curve, the Winter FRF is directly analogous to the Franken FRF. However, the 
frequency scale has been changed since Winter observed a better comparison 
among data when the square root of the diameter was used in the abscissa. 

For internal structure not efficiently coupled to external structure (e.g. , 
internal beams and truss structure) , it is suggested that the acceleration 
spectral density be determined by (a) the solid curve in figure 19 for 
frequencies below 100 Hz, (b) 15% of the solid curve for frequencies above 
500 Hz, and (c) a smooth curve connecting the points at 100 and 500 Hz. 

Winter had added a mass- attenuation factor when heavy equipment is attached 
to the structure. This factor is identical to that described in the Barrett 
method and will be discussed in Section 4.3. 2.2. Winter has also extended 
his method to include vibration induced by aerodynamic noise. This extension, 
suggested by Stewart (ref. 85), is in the form of the efficiency factor shown 
in figure 20. Aer ©dynamically induced vibration is predicted by multiplying 
the acceleration spectral density obtained from the solid curve in figure 19 
(for acoustically induced vibration) by this factor. It is noted that aero- 
dynamic noise is relatively ineffective at low frequencies (usually because 
of a low spatial cross-correlation coefficient) and may be quite effective at 
high frequencies (usually owing to aerodynamic coincidence). 


G”(f) = ACCELERATION 

SPECTRAL DENSITY OF 
STRUCTURAL VIBRATION 

G p (f) = PRESSURE SPECTRAL 

DENSITY OF EXTERNAL 
ACOUSTIC OR 
AERODYNAMIC NOISE 
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RESONANT FREQUENCY 
OF EXTERNAL PANEL 
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f/fl 



Figure 20. - Vibration efficiency of aerodynamic noise relative to acoustic noise. 



4. 3. 1.8 Assessment and Extensional Applications 

The major advantage of the frequency- response methods is the ease in 
calculating the vibration motion of the new vehicle without considering the 
detailed calculations of the structure. Obviously, these methods are highly 
desirable during the early development of a new vehicle, when the structure 
is sufficiently defined or when detailed calculations have not been made. 

The major deficiency of all these methods is their susceptibility to inaccuracy. 
In some cases the overall accuracy is fairly good, while in other cases con- 
siderable overprediction and underprediction have been found (refs. 88 and 
91 through 95). However, inaccuracy should be expected because these methods 
intentionally ignore the very details that determine the response. In 
addition, in those methods which plot their data on a statistical basis, poor 
comparison between predicted and measured vibration is sometimes observed 
because of a too conservative selection of the Statistical parameter. Although 
the various methods have not become standardized on the same unit for express- 
ing the magnitude of random vibration, nor use the same bandwidth, they are 
all convertible to acceleration spectral density, the unit most preferred by 
workers in the field. (E-R-W Method No. 2 is the one exception.) 


Although frequency response methods usually are not used to determine vibra- 
tion stresses, certain modifications to presently available methods can be 
made for their prediction. First, consider a simple mechanical oscillator, 
where the rms stress a ^ in the spring can be related to the rms acceleration 


o.. of the mass: 
w 


a 

s 


(K /4rr 2 f 2 )a 
a n w 


(16) 


where is the ratio of the stress to the (relative) displacement. This 


relationship could be a logical extension of E-R-W Method No. 2, 
2 

values of 6 presented in figure 18. Second, the derivations by 
Ungar (ref. 97), and Crandall (ref. 98) can be used to establish 
relationship between rms velocity and rms stress for distributed 


using the 

Hunt (ref. 96), 
the following 
structures : 


o 
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(EK.K ,/c T ')a. 
<p t h w 
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where E is the modulus of elasticity, 


and K x 

<P 


and , 


the mode- shape factor 


and stress- concentration factor, respectively. For a plate, Hunt shows that 

1.2 < K . < 1.8; from reference 61, 1 < K- < 4. Using K. = %/3 and K- = 1 , 

<p r <t> 1 

Franken compared predicted strain with the strain measured in one- third-octave 

bands on the external structure of the Minuteman missile (ref. 91). The 

strain comparison was only marginal (i.e., measured strain exceeded predicted 

strain by as much as a factor of 5) , even though the acceleration comparison 

was rather good. Further refinements will probably be made as interest grows 

in determining vibration stresses during the early phase of a vehicle- 

development program. 
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4.3.2 Scaling Methods 


When a new vehicle has structural and configurational characteristics and 
fundamental sources of vibration that are similar to the reference vehicle, 
scaling methods make use of experimental data obtained on the reference 
vehicle to predict the vibration on a new vehicle. In general, the two 
vehicles must use the same design philosophy regarding structure, aerodyna- 
mics, and propulsion. Scaling, however, is not restricted to scale modeling, 
which attempts to duplicate all characteristics of the new vehicle. Because 
of the similarities of the vehicles , it would be expected that scaling methods 
are more accurate than frequency-response methods. This is not always the 
case, as is shown later in this section. General scaling methods have been 
used only to predict vibration motion. One reason these methods have not 
been extended to prediction of vibration stresses is that scaling with the 
more accurate scale models have shown only mixed results. These results are 
discussed in the section on model tests (Sec. 5.1.1). 

4.3.2. 1 Condos -Butler Method 

Condos and Butler (C-B) (ref. 99) were the first to utilize data obtained on 
the Titan I missile and the following equation to predict vibration on 
Titan II: 


Sn (£ > ’ I G pn (f)/G pr <f)H<,Sh V < ' 5 h> n ) \r (£) <’ 8 > 

where G^(f) is the pressure spectral density of the acoustic or aerodynamic 

noise averaged over the applicable section of the vehicle; <ph> , the average 
structural mass per unit of exposed area (surface density); and n and r are 
subscripts denoting new and reference vehicles, respectively. Because the 
structures were similar, but were not scale models, C-B did not attempt to 
predict the vibration on a point-to-point basis. Instead, the structures were 
divided into zones of similar location and type of structure. For each zone 
of the reference vehicle, the measured vibration spectra were grouped and a 
statistical analysis made. 

It was observed that the data tended to follow a log-normal distribution, 

which was then used. (This distribution was confirmed later by statistical 

analysis.) To perform the analysis, C-B divided the frequency range of 

interest, 0 to 2 kHz, into 40 contiguous bands of 50 Hz each. In each band, 

the maximum ASD value of each measurement was selected for the analysis, 

rather than the average value. This tended to bias the analysis toward a 

distribution of spectral peaks rather than the distribution of complete 

spectra. Then, using all measurements for the zone, the log-normal mean value 

and standard deviation for each of the 40 bands were calculated, similar to the 

technique used by M-S for octave bands. Finally, the 95%-probability level 

in each band was selected to be used as G.. (f) in equation (18) for predicting 

wr 

the 95%-probability level for the new vehicle. 
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Since the goal of the investigation was the establishment of vibration design 
and test levels for equipment, only measurements on structure adjacent to 
equipment were used* OB selected the following as zones: (1) external 
structure, a stiffened cylindrical shell exposed directly to the pressure 
field; (2) ‘‘outside** structure, secondary structure adjacent to external 
structure but inside the vehicle; and (3) truss structure, a ‘‘space frame** 
consisting of welded beams segments used to attach internal equipment to the 
external stiffeners. The analysis showed no significant differences between 
the three orthogonal directions (supposedly including external structure) , so 
that all directional measurements in a zone were combined. 

To determine the average surface density <ph> , OB used the total mass of 
external and internal structure for the equipment section of the vehicle and 
the exposed area for the section. To use the OB method, it is necessary to 
establish the similarity (although not the identity) of the two vehicles, and 
to determine the required values in the right-hand side of equation (18). 

The zones, type of statistical analysis (or lack thereof), spectral division, 
directions, etc., need not correspond to the parameters used by OB. 

4. 3. 2. 2 Barrett Method 

Barrett (ref. 100) used a method almost identical with that of OB to predict 
acoustically induced vibration on uprated Saturn I and Saturn V vehicles, 
using equation (18) and the vibration measurements from Saturn I. Barrett 
modified the C-B method by using (a) a mass -attenuation factor to reduce the 
predicted vibration when the weight of the equipment was large; (b) a normal 
distribution, modified by a skewness function whose coefficients were exper- 
imentally determined; and (c) a 97.3%-probability level based on the 
calculated distribution. The mass -attenuation factor, to be multiplied by 
the right-hand side of equation (18), is 


Y 


= W / (W + W ) 
n v n e 


(19) 


where and W ^ are weight of the structure and of the equipment of the new 

vehicle, respectively. The 97 .5%-probability level for each zone of the 
reference vehicle was calculated using N zonal measurements and the overall 

th 

rms value X_^ of the spectral envelope of each (or i ) measurement (ref. 101): 


= (50 Hz) 


2050 


2 


G 

max 


< f kh 


( 20 ) 


where G (f ) is the maximum value of the acceleration spectral density 
max x 

2 

(g /Hz) measured in each of the 50-Hz bands (similar to the C-B method) from 
50 to 2050 Hz, and where X. should not be confused with the true standard 

i 
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deviation a.. . = 

wi 
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max 
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L min 


G...(f) df 
wi 
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of the same measurement. Barrett 


calculated the 97. 5%-probability level from 


X = X + no 
cr XX 


where X 


N 

1 $ x. 
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( 21 ) 


(21a) 


(21b) 


(21c) 


(21 d) 


Barrett found that A = 1.64, A = 0.2055, A„ = 0.0155 gave a good experi- 

o 1 z 

mental fit to the 97.5% level. (For a normal distribution where = 0, 

the 97.5% level would be determined at il = 1.96.) Then G.. (f) , to be 

X wr cr 

utilized in equation (18), was determined on the basis that its spectrum 
would have the same general shape as 


N 

WV - s I WVi 


( 22 ) 


but would have a standard deviation equal to X ; that is, 

if 1 

4- J °»r (£) cr d£ (22 ‘> 

f 1 

Figure 21 shows a typical spectrum (f) cr based on this procedure and 
50 vibration measurements made on the aft skirts of the Saturn I fuel tanks. 
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Figure 21. — Vibration-reference spectra for Barrett method on 50 flight measurements on the aft 
skirt of Saturn I fuel tanks. 

Also shown in figure 21 is the mean spectrum G (f ) calculated from 

ITIcLX K. 

equation (22), and the 97.5 percentile G (f ) , determined for each of 

max k. y / • d 

the 50-Hz bands for these 50 measurements, as presented in reference 101. 
Reference 100 shows methods of calculating the average surface density <ph> 
for honeycomb and corrugated structures to be used in equation (18). 


Barrett has also developed a scaling method for mechanically transmitted 
vibration from liquid-fueled rocket engines through structure that is not 
primarily excited by acoustic and/or aerodynamic noise. Engine components 
and structure near engines which do not have large exposed surface areas often 
fall into this category. A scaling equation adapted from Barrett’s work is 


V f > 


N T I W 
n n sp r 
r n 


N T I 
r r sp 



Y G r (f) 


(23) 


where G(f) is the acceleration spectral density of the mechanically trans- 
mitted vibration; N, the number of engines; T and I , the thrust and specific 

sp 

impulse of eacji engine, respectively; and W, the weight of the structure or 
engine component. The value of y is defined in equation (19). Random vibra- 
tion of the rocket engine or of engine-mounted components is however usually 
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influenced by that engine only. Thus, Barrett suggests that = N = 1 be 

used for the combustion chamber, turbopump, and actuator structures and 
components . 

For structure which is excited simultaneously by mechanically transmitted 
vibration and acoustic or aerodynamic noise, such as the structure near the 
engine compartment of the first stage of a launch vehicle, Barrett suggests 
using the following equation to determine the combined effect on the new 
vehicle: 


G .(f) = C G (f) + C G (f) (24) 

nt I wn Z n 

where Ch + C n = 1 , and G.. (f) and G (f) are determined from equations (18) 

1 Z wn n 

and (23) , respectively. Reference values to be substituted into these 
equations, based on measured Saturn I data appear in reference 100. Barrett 
also presents additional reference quantities for sinusoidal vibration, even 
though the measurements showed the vibration to be random. Sinusoidal vibra- 
tion extrapolation to a new vehicle is performed by using the sinusoidal - 
reference quantities of reference 100 and the square root of equations (18) 
and (23), and, if applicable, equation (19) or (24), or both. 

4. 3. 2. 3 Winter Method No, 2 

Winter (ref. 85) has suggested a scaling method similar to his frequency- 
response method discussed previously. In this method, curves for various 
zones are derived similar to those of figure 19, based on vibration and 
acoustic measurements from the reference vehicle. Once these curves are 
determined, the values of surface density and diameter and the predicted 
acoustic environment for the new vehicle are used with these curves to predict 
the vibration for the various zones. 

4. 3.2.4 Assessment 

The major advantage and deficiency of the scaling methods are identical to 
these of the frequency-response methods, even though there are clear structural, 
configurational, and aerodynamic similarities between the reference vehicle 
and the new vehicle. For example, C-B observed as much as a 15-dB difference 
between the predicted and measured 95% levels for the new vehicle (Titan II) 
for the frequency range above 500 Hz, although the difference was usually less 
than 5 dB below 500 Hz. Fortunately, in this case, the prediction was almost 
always greater than the measured percentile. (Since C-B did not use a mass- 
attenuation factor, there was greater overprediction for heavy items of 
equipment.) Similar differences have sometimes been noted for later Saturn 
stages using earlier Saturn measurements and the Barrett method. 



5. EXPERIMENTAL DETERMINATION OF VIBRATION 


Various methods and combination of methods are available to determine design 
adequacy or to predict structural vibration during the various phases of 
vehicle development, or to refine earlier predictions, such as described in 
Section 4. These methods may be divided into the following three categories: 
(1) ground testing, (2) field testing, and (3) flight testing. 

5.1 GROUND TESTING 

Various ground tests may be performed which may have a significant contribu- 
tion to the program. Often the results of these tests are not the sole 
source of data, but must be used in conjunction with the results of other tests 
or analyses, or both. 

5.1.1 Model Tests 


A wide variety of subscale models has been used to establish (1) acoustic 
noise and aerodynamic f luctuating-pressure fields, (2) stability and flutter 
boundaries, (3) propellant-sloshing behavior, (4) natural frequency and modal 
data, (5) structural response to specific dynamic loads, and (6) fatigue and 
other failure mechanisms. The first three uses of models will not be 
discussed in this report except as they may affect the last three uses. The 
primary advantages of using subscale models are the great saving in cost and 
schedule compared to full-scale testing, possible implementation early in the 
program, and often greater accuracy than may be achieved by analysis, 
especially when excessive oversimplification is used. (These physical models 
should not be confused with the mathematical model discussed in Sec. 4.1.1.) 

Models may be divided into the following categories: (a) perfect, or true, 

models; (b) adequate models; (c) distorted models; and (d) dissimilar models. 

5. 1.1,1 Structural-Modeling Theories 

A perfect model reproduces faithfully the significant characteristics of the 
full-scale prototype. Here, the prototype may be the new vehicle or a 
section of the vehicle. However, unless the prototype is quite simple in 
construction (or the scale factor for length is not too large) and no fluids 
are involved, it will probably be extremely difficult or impossible to achieve 
perfect modeling. 

Adequate models give an adequate estimate of the performance of the prototype 
under a given set of conditions without correction, even though certain 
scaling parameters are violated. 

A distorted model requires a correction; but if a correction does not produce 
an adequate estimate, then the model is useless. At present, most structural 
models in use are adequate or distorted. A dissimilar model bears no physical 
resemblance to the prototype, but provides an estimate of performance by 
analogy. An analog computer program set up to represent the prototype is such 
a model. 
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The theory of models requires the selection of certain fundamental units 
(refs. 24 and 102). In structural models, these units are always mass, length, 
and time. If heating effects are important, then temperature is added. Then 
certain independent dimensionless groups tk are obtained by dimensional 

analysis or from mathematical equations which define the dynamics of the 
structure. For structural models, these dimensionless groups are usually one 
or more of the following: the Strouhal, Froude, Mach, and Reynolds numbers, 

or their modification by multiplication with dimensionless ratios. 

The Buckingham pi theorem states that the number of dimensionless groups must 
be equal to the total number of variables (dependent and independent) minus 
the number of fundamental units (refs. 24, 102, and 103). If the problem has 
one dependent variable, the dimensionless group containing the dependent 
variable is a function of the dimensionless groups containing the independent 
variables; i.e., 

17 3 (25) 

Complete similitude is achieved when model and prototype have identical values 

for each ir . . 

x 

There is often more than one dependent variable. In most of these cases, the 
dependent variables are interrelated (e.g., the displacement, acceleration, 
bending moment, and stress at a particular location on the structure), so 
that the remaining dependent variables may be readily determined after the 
first variable is established. In other cases, the dependent variables are 
not related; thus, a more sophisticated approach is required (ref. 104). 

Also, a dependent variable in one problem may be an independent variable in 
another. For example, the natural frequency would be a dependent variable 
when the resonance properties are determined, but would be an independent 
variable when vibration response to a specific loading is obtained. Obviously^, 
the experimenter must be able to distinguish between the dependent and 
independent variables and must assume that the basic phenomenon is independent 
of structural size. 

5. 1.1. 2 Applicability to Space-Vehicle Structures 

For most structural models, scale factors n_^ between the model and the proto- 
type are selected for length, Young 5 s modulus, mass density, and either 

pressure or force; i.e., & = n £ , d =n^E,p =n p, and p = n^p 

r n p £ m p Em p p m 5 p r m 

or F^ = n^ru, respectively. If model and prototype use the same material, 

then n_ = n =1. If it is desired to keep the stress constant, then n^ ~ 1. 

E p ^ p 

The perfect model w r ill have the same mode shape as the prototype in each mode. 

If the model and prototype have their elastic moduli and mass densities 

related by E /p = E /p (i.e., the same speed of longitudinal waves), which 
p p mm 

is often found, then the frequency-scale factor is inversely proportional to 
the length-scale factor n . 

Xj 
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There are difficulties in achieving even a distorted model for most aerospace 
structures subjected to acoustic or aerodynamic noise and/or mechanically 
transmitted vibration. In many cases, it is extremely difficult to scale all 
length parameters adequately. Exceptionally thin skin may be required for the 
model ; thus necessitating special fabrication and handling. The cost of 
scaled rivets and fasteners may be considerable. Since these elements often 
control resonant behavior and determine nonlinear effects, it is important 
that they be installed with great precision. Insufficient room for tools 
often makes this impossible. In addition, surface roughness effects that 
influence friction must also be scaled. 

Problems may arise in producing scaled welded or brazed joints. Composites 
honeycomb, etc., may require special processing to achieve proper scaling. 
Compromises in construction usually have « an appreciable effect on dynamic 
behavior in the higher-order modes, but often little effect in the lower-order 
modes (ref. 102). Unfortunately, for la^ge prototype structures it is the 
higher- order modes that are usually excited by acoustic or aerodynamic noise 
and/or mechanically transmitted vibration. 

If accelerometers are used to measure the vibration of the model, their mass 
loading may have an appreciable effect on the measured vibration magnitude, 
mode shape and natural frequencies, unless they are scaled or their impedance 
small, compared to the driving-point impedance of the structure. The spectral 
density G"(f) for the unloaded structure can however be calculated from the 
spectral density G* (f) measured by the accelerometer (refs. 105 and 106): 

G"(f) = | [Z^f) + Z^<f)]/[Z 1 (f) + Z '2 (f ) | 2 G’(f) (26) 

where Z^ (f ) and Z^(f) are the mechanical (driving point) impedance of the 

structure and the accelerometer at the accelerometer-mounting point, 
respectively; and Z^(f) " For ^ requencies sufficiently below the resonant 

frequency of the accelerometer, Z’(f) ~ i 2irm f, where m is the mass of the 

a a 

accelerometer. 

If propellant sloshing influences the structural vibration in the frequency 
range of interest, another problem is encountered. Slosh natural frequencies 

1 /2 

vary in accordance with (gravity /length) , so that gravity should b« scaled. 

To avoid using an accelerator or centrifuge, most experimenters substitute a 
model fluid which has a greater mass density (ref. 107). Fortunately, the 
frequency range where the effects of acoustic and aerodynamic noise and mechani- 
cally transmitted vibration are predominant is usually above the slosh natural 
frequencies that have an appreciable effect on the structure. However, the 
slosh natural frequencies for both model and prototype must be calculated to 
ascertain this observation. 

When structural models are used to determine the vibration response to acoustic 
noise from rocket or jet engines, best results are obtained when model engines 
are used for generating the acoustic noise. For a perfect model, the 
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fluctuating pressures will be nearly identical to those of the prototype, if 
the model engines use the same propellants and thrust determined from 
2 

T = n n n^>T (refs. 108 and 109). When model engines are not feasible, a ran- 

P Xr p m 

dom acoustic-noise source may be substituted. The characteristics of the 
pressure field may be calculated for the prototype and then scaled to the 
model (ref. 108). 

The vibration response to aerodynamic noise and buffet can be estimated by 
exposing the structural model to aerodynamic flow in a wind tunnel. The Mach 
and Reynolds numbers for the model must be identical with that of the proto- 
type to achieve properly scaled flow conditions. However, no wind tunnel is 
capable of providing simultaneously the desired ranges of Reynolds number and 
Mach number experienced by actual space vehicles. For the forward parts of 
the vehicle, Reynolds -number distortion usually has only a minor effect on 
the spectrum of the fluctuating pressure field (refs. 14 and 15). However, 
little comparative data are available for the aft parts of the vehicle, where 
Reynolds-number differences and the model’s shock reflections off the wind- 
tunnel walls may affect the pressure field appreciably. Unfortunately, most 
wind tunnels have inherently high noise levels, usually caused by general 
turbulence in the tunnel. Also, shocks in the expansion section of a blowdown 
wind tunnel or turbulence generated by the compressor and turning vanes of a 
continuous wind tunnel may add to these high noise levels. If the vehicle is 
aerodynamically c ‘clean,’ 5 this noise can be a substantial contributor to the 
total fluctuating pressure field and the structural vibration of the model, 
thus invalidating the data. If the external shape of the vehicle is 
unfortunately conducive to severe aerodynamic- flow disturbances, such as 
described in Section 3 and shown in figure 3, the aerodynamic noise will 
probably mask the undesirable wind-tunnel noise. 

Occasionally, in wind-tunnel tests, shock waves have been observed to 
oscillate in the transonic regime, as well as to generate increased aerodynamic 
noise. Under certain conditions, the subsonic boundary layer allows the 
higher pressure aft of the shock to feed f orward under the shock. This in 
turn increases the pressure gradient at that point and tends to move the 
boundary-layer separation point upstream. The balance between the initial 
flow conditions that position the shock and the tendency of the shock to move 
forward is unstable. This causes the shock to move randomly about a small 
region of the vehicle so that the flow in this region alternates randomly 
between attached and separated conditions. Oscillating pressure changes as 
high as 0.8 q have been observed in wind-tunnel tests of Apollo and Dyna Soar 
(refs. 15 and 110). A dominant full-scale frequency range of 0 to 2 Hz has 
been observed. However, shock oscillations have not been observed during 
flight. The increasing speed of the vehicles does not permit the shocks to 
c ‘stabilize 5 5 in the transonic regime, and causes the shocks to pass through 
those conditions conducive to shock oscillations. Thus, oscillating shocks 
could occur for a vehicle operating at constant speed in the transonic regime, 
such as in a wind-tunnel test, but would not for a vertically rising vehicle 
following a conventional trajectory. 

Structural buckling of the model under static loading may be another problem 
if the prototype structure normally stands vertically and the model is tested 
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horizontally in the wind tunnel. Experimental data on the buckling of 
structures similar to the model may be used to determine if horizontal 
testing will be a problem (refs. Ill through 113). (Results of theoretical 
analyses on idealized structures should be used with great care because they 
often provide an un conservative estimate of critical buckling loads.) 

Mechanically transmitted vibration is not often modeled. Best results can 
probably be expected using model sources for the excitation. If the source 
of mechanically transmitted vibration is the rocket engine, it will probably 
be difficult or expensive to obtain a model rocket engine with sufficient 
detail to provide proper dynamic-force scaling. Other sources, such as pumps 
and generators, can be expected to have similar difficulties. 


If an electrodynamic or hydraulic shaker is substituted for the model source, 
it should apply a spectral density G^f^) , which can be calculated from 

equation (26) after certain vibration and mechanical impedance measurements 
are made on the prototype source and on the adjoining test structure (e.g., 
the full-scale rocket engine and its test stand), and on the model structure 
(ref. 105). The acceleration spectral density G^(f^) to be used in the 

right-hand side of equation (26) can be calculated from the acceleration 
spectral density G* (f ), which is measured at a location near what is normally 


the interface of the prototype source with the prototype structure, t during 
excitation of the source (e.g., static firing of the rocket engine): 

1 3 


G'(f ) = (n.niti|[/DJG»(£ ) 
mm £ E p' p p p 


(27) 


The frequency should be scaled in accordance with 

2 . 

f = n (n /n ) 2 f 
m Z p/ E p 


(28) 


The mechanical impedance Z^(f^) and vs? 9 can ca l cu ^ atec * f rom the 

driving-point impedance measurements Z.j (f^) and Z^Cip), made separately on 

the prototype source structure and the test structure, respectively, at the 
interface: 


2 i» ( y - vy/^vi' 
z vy - ^y/^v* 


(29) 

(29a) 


The mechanical driving-point impedance Z/J^f^) can be determined by direct 

measurement of the model structure at the normal interface with the shaker 
which represents the model source. Impedance-measurement techniques are 
discussed in references 114 through 118. 

The vibration response to acoustic or aerodynamic noise, or to mechanical 
transmission from the source, can be scaled correctly only if the structural 
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damping is not affected by the scaling. As discussed in Section 4. 1.2. 2, 
damping occurs in one or more of the following three forms: material damping, 

friction damping, and acoustic radiation. Material damping is dependent upon 
the mode shape and the maximum vibration stress in each mode. Although 
figure 6 shows that material damping is nonlinear, the material damping of 
the model will be identical to that of the prototype if (a) the mode shape is 
identical, (b) the model and prototype structures are made of the same 

material, and (c) ru, = 1 . If the maximum stress is less than the endurance 

P 

limit for both model and prototype structures, the resonant magnification 
resulting from material damping changes only slightly for different stress 

levels, obviating the need for n_ = 1 . 

P 

An analysis performed by Greenspon, however, has shown that it is almost 

impossible to scale material damping perfectly (ref. 119). He assumed thermal 

relaxation to relate the internal friction to the temperature rise of material, 

and found that the thermal conductivity k and the specific heat c of the 

s s 

material must be scaled in accordance with 

"k = VWX (30) 


When the same material is used for model and prototype (n = n = n = n 

K.S C S ill p 

= 1), equation (30) is violated unless n^ = 1 (i.e., no length scaling). Even 


when different materials are used, in most cases they have almost the same 


speed of longitudinal waves (i.e., n 


V' 


Thus, under these conditions, it 


would be fortuitous to find a model material that would satisfy equation (30) . 
It should be noted that the scaling difficulty occurs in the heat conduction. 


After the heat is conducted to the surface, it is transferred to the surround- 
ing fluid by conduction and convection, and to the other surfaces and spaces 
by thermal radiation. For scaling the heat conduction across the layer of 
fluid adjacent to the structure, the model fluid must be selected to provide 
the same Reynolds number as the prototype fluid; that is. 


(p_ w l /\x ) = (p_ w £ /y ) 
fp p p p fm mm m 


(31) 


where p^ and y are the mass density and viscosity of the fluid, respectively, 

and w is the oscillatory velocity of the fluid at a point on the structural 
interface. Even so, the thermal conductivity and specific heat c^ (at 

constant pressure) must also be scaled in accordance with equation (30) . If 
the same fluid (usually air) is used for model and prototype, then 

=n = n = n = 1 , again violating equation (30) unless n = 1. It 

Cf E p ** 

would be again fortuitous to find a model fluid that would satisfy equation 
(30). Thermal radiation will be scaled if all spatial relationships are 
scaled in accordance with n^ and the following characteristics of the model 

structure and all adjacent surfaces are identical to those of the prototype: 
materials, temperature differences, emit.tances, and absorptances (ref. 120). 
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This scaling would be difficult to accomplish practically, even though the 
experimental data of figure 6 does not shew the scaling to be as great a 
problem as is indicated by the theory. 

Because friction damping is controlled by the kinetic coefficient of friction 
between the mated surfaces, it would appear on a theoretical basis that scal- 
ing should be no problem despite the fact that it is a predominant source of 
nonlinear vibration. Since the static friction must be overcome for slip or 
sliding to occur, (a) the static as well as the dynamic coefficient of 
friction for the model and prototype structures must be identical, and (b) the 

fluctuating pressure fields must be scaled in accordance with n^ = 1 (refs. 47 

P 

and 49). However, the first condition is usually impossible to achieve in 

practice. Acoustic radiation from the model and prototype will be identical 

if the spatial conditions and adjacent structures are scaled and if 

(c /cl ) = (c / c ' ) and (p £ /p ) = (p_ /p ), where c is the speed of sound 
ap Lip am Lm rp p tm m a 

in the fluids; c' , the speed of longitudinal waves in the material; and p and 
L 

p^, the mass density of the structure and fluid, respectively. If the same 

material and fluid are used for the model and prototype, proper scaling should 
be achieved. 


For viscoelastic damping, it can be shown that model and prototype have 

identical resonant magnifications if the model uses the same materials (metal 

and viscoelastic) as the prototype. To achieve proper scaling for air 

pumping, it is necessary that (c /cl ) = (c /cJ ) and (P /P ) = (E /E ) = 

ap Lp am Lm op om pm 

1/n^, where P q is the absolute ambient pressure. Unfortunately, it would be 

extremely difficult to satisfy the latter relatipnship , except under special 
circumstances, and impossible to satisfy if the same materials and fluid were 
used for the model and prototype. However, theoretically-determined 
corrections can be applied to model results to predict prototype performance 
accurately . 


Scale models can also be used to predict fatigue behavior. Obviously the 

behavior is dependent upon the material, so that the model and prototype must 

be made of the same material (i.e., n_ = n =1). Since the number of cycles 

L p 

to failure is dependent on the stress, the scale factor n^, = 1 must be used. 

If scaling is perfect, the time- to-f ailure T^ for model and prototype will be 

related by T^ = n^T^. However, fatigue tests performed by Gray on simple 

panels exposed to jet noise, using models of various scale factors n^ , show 

imperfect scaling in the time to failure (ref. 121). From these tests, Gray 

X 0 2 

observed by T- = n_T_ , where X = 1/n ' , indicating that his models produced 
fp 36 fm 36 

conservative estimates for prototype-f atigue life. In most cases, this result 
can be explained by lack of scaling of the nonhomogeneity of the material 
(grain size, microscopic flaws, stress concentrations , etc.) between model and 
prototype. The effects of size on fatigue is discussed in reference 122. 
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Scale models are seldom used for other vibration failure mechanisms; for 
example, collision between adjacent structures, wear, and operational drift. 

5. 1.1. 3 Assessment 

To be of practical utility during the development of a space vehicle, a 
subscale structural model must provide the prediction of prototype performance 
with sufficient accuracy, with or without correction. A review of published 
literature on the comparison of model and prototype data shows that whc n the 
prototype was a highly complex structure, the model did not provide an 
adequate estimate of prototype performance, but when the prototype was fairly 
simple, the model did provide an adequate estimate. For example, acoustic 
noise was applied to various external panel sections of the highly complex 
Snark missile and two quarter-scale models (i.e., n^ = 4) , and the vibration 

measured at many locations throughout the three structures (ref. 123.). The 
original model had many deficiencies in structural details, although most of 
the basic features of the prototype were properly scaled. The improved model 
was an attempt to correct most of the important deficiencies. Figure 22 shows 
a typical comparison of one- third-octave vibroacous.tic- transfer functions 
(ratio of rms -vibration acceleration to rms-acoustic force) versus scaled 
frequency for these structures. 
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Figure 22. - Comparison of vibroacoustic-transfer functions for Snark missile structure and two 
one-quarter scale models. 
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The prediction produced by the original model was almost always in error, 
often by more than 20 dB. The improved model provided some improvement in 
accuracy, as shown in this typical figure; but the error was still too large 
(occasionally in excess of 20 dB and often in excess of 10 dB) to preclude 
its practical utility or to justify its cost. Another example of inaccuracies 
was observed between the Apollo boilerplate service module and its one-tenth 
scale model, shown in figure 23 (ref. 124). Great pains were taken to scale 
nearly all structural details in the model. Transonic flight data for the 
full-scale boilerplate were compared with wind-tunnel data for the model. 
Figure 24 shows a reasonably similar shape between the two normalized random 
strain spectra. Unfortunately, the model underpredicted the rms value of the 
strain by 8 dB during the most severe transonic period. On the other hand, 
the model underpredicted the rms acceleration of the flight prototype by 
only 2 dB. Unfortunately, the model missed the predominant frequency range by 
a factor of 2, as shown in figure 25. Reasons for the inaccurate predictions 
of the Snark and Apollo models are discussed in references 123 and 124. 

Similar results were noted on a Nimbus spacecraft at higher frequencies (ref. 
125). 

On the favorable side is the comparison of results for simple panels exposed 
to jet noise, mentioned previously in reference 121. Figure 26 shows the 
excellent prediction of the prototype-strain spectral density from the 
one-sixth-scale model data. Similarly, models of simple shells (cylindrical, 
conical, and hemispherical) have shown excellent estimates of prototype 
natural frequencies and mode shapes for lower-order modes. However, no effort 
has been made to determine the degree of complexity obtainable before 
prediction inaccuracies occur. 

5.1.2. Vibration Tests 


5. 1.2.1 Test Classification and Philosophy 

Laboratory vibration tests may be categorized as follows: (a) design develop- 

ment, (b) qualification, (c) acceptance, and (d) others. Design-development 
tests are performed to provide certain information to the designer regarding 
the modal characteristics of the structure (ref. 33) , the vibration response 
to certain applied loadings, the adequacy of the design, and, in some cases, 
the vibration-failure mechanisms, locations, and levels. Qualification tests 
are performed on flight -quality hardware to demonstrate the adequacy of the 
design and fabrication methods for flight, usually before the first flight. 
Acceptance tests are performed on articles intended for flight usage to 
demonstrate that a predetermined minimum adequacy had been achieved during 
fabrication and assembly. Any of these may be performed as vibration tests 
combined with other loads and environments. Other tests include (a) reliability 
tests, to demonstrate the variation of the failure mechanism between items of 
hardware under prescribed types of loading; (b) fragility tests, to map the 
failure threshold perimeter under a variety of loadings or frequencies; and 
(c) certain interim tests between development and qualification. There is no 
common terminology used throughout the aerospace field to describe a particu- 
lar type of test, so it is not uncommon in two different vehicle-development 
programs or in two different organizations to have two names for the same 
type of test. 
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STRAIN SPECTRAL DENSITY 





Figure 26. - Comparison of vibration-strain spectra for prototype and model panels excited by jet noise. 

Usually, the selection of design-development tests and their test conditions 
are at the option of the designer and/or program management. Qualification 
tests are nearly always required by contract, with test conditions usually 
specified as somewhat more severe than those anticipated for flight. Quali- 
fication tests are often waived when the structures are too large to be 
tested. In these cases, the structures are often { * qualified 9 9 by vibration 
analyses and/or tests on certain structural sections. Laboratory vibration- 
acceptance tests may or may not be required by the contract. If required, 
the test conditions must be selected with great care, since the test degrades 
the hardware by exposure before flight, while the flight adequacy is being 
demonstrated simultaneously. Hardware used for testing, other than accept- 
ance testing, is never used for flight, except under unusually special 
circumstances. 

Laboratory vibration tests are nearly always performed on equipment, sometimes 
on secondary structure (almost always with equipment items or dummies), and 
occasionally on primary structure. 

The subject of equipment testing will be limited in this report to the 
discussion of design and test criteria. However, it is likely in the future 
that some equipment items may serve simultaneously as minor or major load- 
carrying members, so that this division may eventually disappear. 

There is no uniform philosophy throughout the aerospace industry regarding 
the types of primary or secondary structures to be tested. Most spacecraft 
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contracts require vibration testing of the entire payload. For example, 
figure 27 shows the Ranger spacecraft excited by an electrodynamic shaker 
simulating the longitudinal vibration applied by the launch vehicle during 
the launch phase of the mission (ref. 126); figure 28 shows a similar test 
on the Gemini spacecraft (ref. 127). Figure 29 shows the Surveyor spacecraft 
excited by three electrodynamic shakers simulating the vibration applied by 
the three vernier rocket engines during the lunar- descent phase of the mission 
(ref. 128). The spacecraft is inverted to permit operation of the radar sub- 
system in the free field. 

Only a few contracts require vibration testing of launch vehicles or large 
spacecraft because most are too large and too heavy to be tested conveniently 
or economically. In some cases, sections of the launch vehicle are tested; 
figure 30 shows the vibration- test setup for the thrust structure and aft 
skirt of the Saturn S-II stage (ref. 129). Eight hydraulic shakers excited 
the structure through a large fixture. Dummy engines and equipment, simula- 
ting the weight and center of gravity of the actual items, were used to 
provide a simulated loading for the primary and secondary structure of the 
section. 


5. 1.2. 2 Test Facilities 

There is a wide variety of test equipment available for performing laboratory 
vibration tests. The electrodynamic shaker is the most popular because of 
the ease of controlling the test. It is used mostly in the 5- to 2000-Hz 
frequency range, although the useful upper- frequency limit usually depends 
upon its force capacity and size. Electrodynamic shakers with force- genera- 
ting capacities up to 30 000- lb rms are commercially available. 

Hydraulic shakers are used mostly in the 0- to 500-Hz frequency range, although 
the useful upper-frequency limit is dependent upon the design features of the 
shaker, which vary widely, as well as on its force capacity and size. Hydraulic 
shakers with force-generating capacities up to 250 000-lb rms are commercially 
available. 

Test systems for providing power and control to these shakers also vary widely. 
For performing sinusoidal vibration tests, a typical system has an electronic 
oscillator driving an electronic power amplifier through a gain-control cir- 
cuit, as shown in figure 31. The power amplifier in turn drives the shaker 
(actually a translational motor) which is attached to the test fixture and 
test specimen. 

In nearly all tests, a control accelerometer is attached to the test fixture 
at a point adjacent to the mounting point of the specimen. If there are two 
or more mounting points, the accelerometer is located adjacent to one of them. 
The accelerometer signal is either fed directly to the gain control circuit, 
so that the shaker is automatic-gain-controlled (AGC) , or fed to a meter, 
which is read manually, and the gain- control circuit is then adjusted 
manually. 
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Figure 27 r Laboratory vibration test of the Ranger spacecraft. 
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Figure 30 -Laboratory vibration-test facility of the thrust structure and aft skirt of the S-ll stage, 
Saturn V launch vehicle. 
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TEST SPECIMEN 



Figure 3 Vr Simplified schematic diagram of a typical system for performing sinusoidal-vibration tests. 

To perforin random-vibration tests, a typical system contains a random-noise 
generator driving an electronic power amplifier through a spectrum equalizer, 
which consists of a large number (usually 40 or 80) of contiguous bandpass 
filters, as shown in figure 32. Each filter covers a narrow band portion of 
the spectrum (usually 50 or 25 Hz) . The gain of each filter is varied auto- 
matically by means of AGC, or manually. The output signals from all the 
filters are summed and fed to the power amplifier, which in turn drives the 
shaker, fixture, and specimen. The control accelerometer is attached to the 
fixture at a point adjacent to the specimen, as described previously. The 
accelerometer signal is fed to a spectrum analyzer, which consists of a set 
of filters that have bandpass frequencies identical to those in the spectrum 
equalizer. The signal from each filter in the analyzer is either fed to its 
respective equalizer circuit, controlling the equalizer gain by AGC, or fed 
to a meter (one of 40 or 80), which is read manually. The gain of all 
equalizer circuits is adjusted manually. These systems are obviously expen- 
sive. It must be noted that figures 31 and 32 show only the critical c ‘black 
boxes 5 5 of the test equipment. In actuality, there are many additional items 
of equipment needed for control and overload protection. 

5. 1.2. 3 Qualification-Test Requirements and Their Selection 

Because of its contractual status, the qualification test is probably the 
most important type of vibration test performed during vehicle development. 
The intent of most qualification tests is to accept flight -adequate hardware 
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Figure 32 “Simplified schematic diagram of a typical system for performing random vibration tests . 
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and to reject flight -inadequate hardware. Short of actual flight test, it is 
impossible to achieve this goal in practice. Thus, it is necessary to consider 
the tradeoffs between undertesting and overtesting. Obviously, it is easy to 
avoid nearly all undertesting by increasing test levels until only the most 
rugged hardware survives. In this case, much of the hardware that is really 
adequate for flight must be redesigned and retested, often causing unnecessary 
weight increases and cost and schedule slippage. Many test specifications are 
written by personnel whose main assignment is to ensure high reliability and 
structural integrity, leaving to program management and other groups the weight, 
cost, and schedule problems resulting from overly conservative test levels 
(ref. 130). With this information, it is not surprising to find that there 
are many instances where program management has decided to fly space vehicles 
with hardware that has failed vibration-qualification tests. 

The vibration problem that has probably been discussed in the aerospace 
industry more than any other over the years is that of the selection of quali- 
fication test levels and durations. The factors that should be considered 
in this selection include: 

• The degree of similarity between the test and the flight 
configuration. 

• The characteristics of the test facility. 

• The method of controlling the test. 

• The degree of confidence in the flight vibration 
predictions . 

• The number of test specimens to be used in the 
qualification of each hardware item. 

• The economic balance between the cost of possible flight 
failures and the cost of possible test failures. 

The degree of similarity between test and flight configurations will have a 
major influence on the lowest frequency for which valid testing can be per- 
formed. Often, this frequency limit can be readily determined by comparing 
the modal characteristics at the lower resonance of the test specimen with 
the modal characteristics expected for the flight configuration. If these 
occur at significantly different frequencies, or if no similar mode shapes 
exist, no valid vibration testing is possible for these modes (ref. 131). 

For this reason, it is common practice to test a spacecraft by ‘ ‘notching’’ 
the vibration test level in the vicinity of its fundamental resonance. 

Actually, an identical natural frequency or mode shape could exist only by 
coincidence because of the different mechanical impedance characteristics 
between the test fixture and the flight structure adjoining the hardware 
section under consideration. It is possible, however, to force the specimen 
to high levels at frequencies corresponding to the flight-resonant frequencies 
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and c ‘notch 5 * the levels at frequencies corresponding to the test-resonant 
frequencies and flight-antiresonant frequencies, if the shaker capacity is 
large enough and if the resonant characteristics of the flight configuration 
(usually not yet flown) are known with sufficient accuracy in the frequency 
range of interest. However, this is almost never done. 

For smaller structures, such as the spacecraft shown in figure 27, it is 
common practice to vibrate the specimen by a single shaker through a test 
fixture. There are usually great pains taken to restrain the vibration motion 
of the shaker and fixture to a single direction, and to make the test fixture 
as rigid and as massive as the force limitations of the shaker permit. Such 
care helps avoid any contribution to the resonant behavior of the specimen by 
the fixture, except at higher frequencies. This is to be compared with 
the flight configuration where the adjoining structure is part of the resonant 
behavior, often allowing motion in all three orthogonal (mutually perpendicu- 
lar) directions. This permits multidirectional resonance, and, in the case 
of most hardware that is attached to adjoining structure through two or more 
points of attachment, allows variation of amplitude and phase between attach- 
ment points throughout the frequency range. Thus, there is a considerable 
lack of simulation in using single- shaker restrained-motion vibration testing 
with rigid fixtures. Much of this problem can be reduced by using part of the 
flight structure adjoining the test specimen as a portion of the test fixture. 

The variation between attachment points can often be improved by multiple 
shakers under separate control, but this increases the cost of testing. In 
addition, the shaker applies the vibration force in only one direction, whereas 
the flight excitation is applied in the three orthogonal directions simul- 
taneously. Usually, it is assumed that vibration applied sequentially in each 
of three orthogonal directions for a specified time duration per direction has 
the same damage potential as vibration applied in al ; three directions simul- 
taneously for the same duration. However, this assumption can only be valid 
for resonant behavior in one direction and by coincidence for multidirectional 
resonance . 

As stated previously, most vibration tests are controlled from a signal 
generated by the control accelerometer, which is attached to the test fixture 
adjacent to the mounting point or to one of the mounting points. In the large 
majority of tests, the specimen interfaces with the fixture at several points 
or along lines or surfaces. Multidirectional resonant behavior of the test 
specimen (including rotation) , which is coupled with fixture and shaker 
resonances at the higher frequencies, usually causes the vibration amplitude 
to vary between mounting points or along the interface. Since the control 
accelerometer senses motion at only one of these points, undertesting or 
overtesting usually occurs at the other mounting points, the magnitude of 
which depends on the mode shapes of the multidirectional resonances. To avoid 
this undertesting or overtesting during a sinusoidal vibration test, several 
control accelerometers (usually from three to six) can be used and their 
signals monitored simultaneously to select the lowest or highest amplitude 
for test control. Special test equipment is commercially available for imple- 
menting this selection automatically. 


73 


Another technique employs the averaging of the signals from the several 
accelerometers. Signal averaging may be used for sine or random vibration 
testing. Signal averaging may be performed in one of two ways: (1) instan- 

taneous averaging and (2) time division multiplexing (TDM). Overtesting can 
occur with the instantaneous averaging technique in any mode when two of the. 

control accelerometers are mounted at locations that are vibrating 180° out 
of phase. This problem is especially acute when only two accelerometers are 
used, or when the out-of -phase accelerometers are vibrating at high amplitude 
relative to the other accelerometers. This problem may be avoided by using 
TDM, which samples the accelerometer signals in sequence to form a composite 
signal which is then used for manual or automatic test control. However, 

TDM is not without problems, such as the selection of sample time durations 
from one signal to the next, although these effects can be minimized (ref. 
132). Special equipment is commercially available for implementing test 
control with TDM. 

As seen from the time histories and the changes in the vibration spectra 
illustrated in figure 1, the flight vibration is nonstationary, especially 
for the period from liftoff through q raax « However, most presently available 

vibration test systems are capable of controlling stationary tests only. 

Thus, stationary testing is employed for nonstationary flight conditions. 
Several test techniques may be used, but all permit overtesting in the time 
and/or frequency domain. - The most popular is the use of a vibration spectrum 
which envelopes the vibration spectra for the various prelaunch and mission 
events that exhibit high vibration (acceptance vibration testing, static 
firing, liftoff, transonic, Q 7nax ? stage-powered flight for hardware near 

engines, etc.). Often the time duration for this stationary laboratory test 
is selected to be the sum of the € c ef f ective* 5 durations of these prelaunch 
and mission events. Thus, considerable overtesting may result compared to 
flight conditions. Overtesting can sometimes be minimized when qualification 
testing is performed by using several tests in sequence, each test represent- 
ing a different phase of the mission or of prelaunch condition, a.nd having a 
different vibration spectrum. However, this series of test setups, checkouts 
etc. , may be tedious to perform. In the future, special equipment or tech- 
niques will probably be available to permit nonstationary vibration testing. 
Two seem most likely: (1) the use of a nonstationary signal, which has been 

prerecorded onto magnetic tape that is fed directly to the kind of electronic 
power amplifier shown in figure 32 (using no AGC) , after the system has been 
calibrated manually with a low-level white random spectrum; and (2) the use 
of an electronic programmer between the control accelerometer and AGC 
(fig. 32). 

The degree of confidence in the flight-vibration prediction should also 
influence the selection of the qualification test levels and durations. 
Classical techniques of vibration prediction applied to a sophisticated 
representation of the flight structure and test specimen, discussed in 
Section 4.1, should permit a rather accurate selection of test levels up to 
a frequency where the representation becomes increasingly less accurate. At 
higher frequencies the test levels must be increased, usually in some arbi- 
trary manner to avoid the possibility of undertesting. Similarly, the 
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statistical- energy analysis (SEA), discussed in Section 4.2, may be used for 
selecting qualification- test levels at higher frequencies. However, since 
SEA is less precise at high frequencies than classical, techniques are at low 
frequencies, the test levels derived by SEA must be somewhat increased to 
avoid undertesting. The amount of the increase will also depend upon the 
degree of sophistication used in the structural representation. If sufficient 
detail is included, the test level should not be more than the SEA upper 
bound described previously, which was established by setting the structural 
damping to zero and utilizing only the coupling loss factor. 

Extrapolation techniques (Sec. 4.3) are most often used for test level selec- 
tion. The degree of confidence in the prediction will depend upon the 
similarity between the reference and the new vehicle, and how the reference 
vehicle data were grouped, and if applicable, statistically analyzed. If 
the reference and new vehicle structures are quite similar, it might be 
assumed that the vibration prediction would be quite accurate. However, as 
discussed previously in Section 4.3.2 on scaling methods. Condos and Butler 
found considerable differences between vibration predictions and subsequent 
measurements on a new vehicle, even when excellent similarity existed between 
the reference and the new vehicle, and the reference vehicle data were grouped 
into zones for similar location and type of structure. From this, as well as 
from similar experience on the uprated Saturn I launch vehicle, increased 
test levels to avoid the possibility of undertesting seem in order. 

It would appear that a greater increase would be in order for frequency- 
response methods (Sec. 4.3.1), since even less similarity between vehicles 
is generally found. However, in all extrapolation methods, the number of 
measurements, the grouping of the data, and the type of statistical analysis 
(if used) often has a large influence on the vibration level predicted and 
thus on the resulting test level. In regard to the number of measurements, 
for example, Stevens compared eight vibration measurements made on 
Apollo- service-module external panels (considered simple aerospace structure) 
with the total population of 180 measurements made on these panels in the 
radial direction when the service module was acoustically tested (ref. 133). 
These eight measurements were selected because their locations corresponded 
with flight-measurement locations. Figure 33 shows that generally there was 
a 3-dB difference between random vibration spectral envelopes for the small 
and large number of measurements. This comparison has a direct bearing on the 
selection of test levels when only a small number of measurements are used 
for a zone of the reference vehicle in a particular direction. 

In regard to data grouping, care is usually required (ref. 86). For example, 
for statistical analysis, it is common practice to put all vibration measure- 
ments for a particular zone of the reference vehicle into the same group, 
independent of the direction of the measurement. For structure exhibiting 
vibration predominantly in one direction, this type of data grouping can 
produce two effects on the test level. First, overtesting can occur in the 
two directions of lower vibration when the test level is selected on the basis 
of the entire population of three directional measurements. Second, over- 
testing can even occur in the direction of highest vibration if the test 
levels are selected on the basis of a statistical analysis which uses the 
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Figure 33.- Comparison of spectral envelope for 8 and 180 vibration measurements made on the 
external panels of the Apollo service module during acoustic testing . 

mean value and standard deviation of the entire data population. A few data 
points for the two directions of lower vibration can cause a significant 
increase in the standard deviation without an equivalent decrease in the mean 
value. Thus, if the test level was based on a certain high probability level, 
such as the 95 percentile, a higher test level would be produced by using all 
measurements (independent of direction) than would be produced by using only 
measurements for the direction of highest vibration. 

Although not often considered, the number of test specimens used for qualifi- 
cation testing should influence the selection of the test level. Of particular 
concern is the variation of the fragility level from specimen to specimen. 
Fortunately, the primary failure mechanism of primary and secondary structure 
is fatigue caused by response at the lower resonant frequencies. The fragility 
(time- to-failure) variation is often acceptable here because manufacturing 
tolerances can be reasonably controlled. The opposite is sometimes observed 
for equipment items, since the predominant fragility range can occur at the 
higher frequencies where manufacturing variations cannot be adequately 
controlled. This is because tolerances are usually controlled to mils (i.e., 
milli-inches) , whereas relative vibration displacements often occur in the 
micro-inch range, and also occur where other failure mechanisms than fatigue, 
such as threshold crossings, may be dominant. When test specimens are 
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Figure 34r Predicted vibration spectrum and a five-segment spectral envelope. 

antiresonance (refs, 136 and 137). Thus, tolerances are usually specified 
for the test spectrum, often in the vicinity of ±3 dB for narrow-band 
analysis. If it is desired to avoid any possibility of undertesting, the 
test spectrum must be increased to match the lower tolerance. Then additional 
over testing is produced at other frequencies. However, there has been much 
discussion regarding the necessity for this increase, since many test failures, 
as well as flight failures, are likely to occur at frequencies corresponding 
to antiresonances, as measured by an accelerometer on the test fixture or at 
the flight interface, although these antiresonances would not occur at 
identical frequencies (refs. 138 and 139). 

Many of these problems may be avoided by using a test method developed by 
Curtis (refs. 140 and 141), sometimes called 5 ‘random on random*’. Three 
high-level random vibration narrow-band ‘‘spikes 55 are applied simultaneously 
with a lower-level wideband random vibration spectrum, as shown in figure 35. 
Curtis suggests that (a) the levels for the wideband and narrowband portions 
be determined from the mean and 98% values, respectively, which are calculated 
from equation (15) in Section 4.3.1; (b) 12.5% bandwidths be used for the 
three narrow-band spectra; and (c) the ‘ ‘spikes 5 5 be stepped sequentially over 
a given frequency range (fig. 35), but never overlap. Of course, other 
variations are possible. 

Undertesting can result when a possible failure is dependent upon simultaneous 
high-level excitation at two or more frequencies, which might occur in flight 
but is not programmed into the narrow-band sequence. Overtesting and/or 
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Figure 35. - Random vibration-test spectrum consisting of high-level narrow band "spikes" 
superimposed on lower-level wideband spectrum. 

undertesting can also occur, depending upon the time duration of application 
of each narrow-band combination, as well as the wideband spectrum. However, 
Curtis reports good correlation between flight and laboratory results for 
equipment used in aircraft. 

5. 1.2. 4 Alternatives for Insufficient -Force Capacity 

If the test level is too high or the test specimen too heavy, the shaker’s 
force capacity may be exceeded. The obvious solution is to use a larger 
shaker or multiple shakers. But in many cases, the largest commercially 
available shaker is already considered (for the frequency range desired) or 
spatial limitations preclude the use of a sufficient number of shakers. In 
these cases, alternate measures are required. First, the Curtis method should 
be considered, since a much lower rms acceleration is required, compared to 
that for the overall envelope of the high-level narrow-band 4 ‘spikes.’ 5 
Second, consideration should be given to using high-level narrow-band random 
vibration which is applied simultaneously at the predominant resonant 
frequencies of the test specimen for the desired time duration, with the sur- 
plus shaker capacity (if any) being used to produce a lower-level wideband 
random spectrum outside the resonant frequency ranges. A low-level sinusoidal 
resonance search may be used to establish the predominant resonant frequencies. 
Third, a narrow-band random spectrum swept across the frequency range may be 
used, as described by Booth and Broch (ref. 142). Finally, a sinusoid may be 
used for resonant dwell at the various resonant frequencies, or it may be 
swept across the frequency range. Some type of sine-random equivalence is 
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usually assumed (e.g., ref. 66), although it is unlikely that it can be 
achieved (ref. 143). Of these substitutes, the Curtis method is probably the 
most realistic; the random-at -resonance method, the most conservative; and 
swept random, swept sine, and sine dwell, the lowest in cost. 

A completely different technique is to specify and control force, rather than 
control acceleration in vibration testing. Although this form of testing is 
harder to implement, and vibration force measurements are not normally made 
in flight, it is possible that much more realism may be achieved with force- 
control testing, thus avoiding under testing and over testing. Otts and his 
associates have demonstrated that random-force vibration testing can be 
implemented practically in the laboratory (refs. 144 and 145). 

5. 1.2. 5 Acceptance -Test Requirements and Their Selection 

Acceptance testing of articles to be used later for flight is seldom performed 
on structure, but is often performed on equipment items. Variations in 
structural integrity are usually small for structures, especially in the 
lower-order modes, so that qualification testing of an earlier specimen at 
higher levels (to avoid the possibility of undertesting) , plus normal inspec- 
tion of the flight-usage structure, should assure the adequacy of the 
structure for flight. On the other hand, there can be large variations in 
the structural integrity and operational performance of equipment, usually 
because the predominant response is in the higher-order modes. Thus, quali- 
fication testing at higher levels, plus normal inspection, will not necessarily 
assure flight adequacy. As stated previously, the test levels and durations 
must be selected with great care so that degradation resulting from preflight 
exposure is minimized while flight adequacy is being demonstrated. 

Most of the problems found in the selection of qualif ication-test levels are 
also encountered in the selection of acceptance- test levels. However, the 
6 * solutions 5 9 to these problems are usually different. For qualification, a 
series of test inadequacies are usually countered by a series of test-level 
increases (i.e. , overtesting) , whereas for acceptance, the test inadequacies 
are usually judged equally for their effects on undertesting and overtesting. 

It should be remembered that the cost of an acceptance-test failure is 
almost always much less than the cost of a flight failure. It is not un- 
common to find acceptance-test levels and durations specified as a certain 
percentage of qualif ication-test levels and durations. In some cases, stage 
static-firing tests (i.e., preflight acceptance firings of rocket engines 
in their stages) are used in lieu of laboratory acceptance tests. This is 
usually justified by assuming that there are fairly minor differences between 
flight and static-firing vibration levels, but this assumption is often 
violated. For example, for certain sections of the upper stage of a launch 
vehicle, the acoustic noise and the resulting equipment vibration from the 
upper-stage rocket engine (s) during static firing may be much less than the 
noise and vibration from the first-stage rocket engine (s) during liftoff. 

Consideration must be given additionally to the differences in the cost 
between laboratory acceptance- test failures and static-firing failures. It 
is thus easy to see why acceptance tests are not always used and why there 
is no uniformity in acceptance- test level selection. 
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5. 1.2. 6 Design- Development Test Selection 

There is a tremendous variety of design-development tests (DDT) , since the 
major goal of these tests is usually to provide information to the designer 
and program management on the. adequacy of the design and on the vibration 
characteristics of the hardware, in case the design is inadequate. These 
tests are usually performed quite early in the vehicle development so that 
program schedules are not affected if design changes are considered necessary. 
The test-to-f ailure technique is often used during design development to 
establish margins over flight and qualification- test conditions. Since they 
are performed to acquire design information, there is no stigma attached to 
a DDT failure. It is not uncommon to find that the first items off the 
production line are used for vibration DDT, since they often contain minor 
flaws, patches, etc., which make them unsuitable for other uses, but they 
are often still adequate for DDT purposes. Program management should consider 
DDT as an inexpensive tool in avoiding costly failures later in the program. 

5. 1.2. 7 Assessment 

It can be concluded from the preceding that it is impossible to perform a 
vibration test which, even under the best circumstances, avoids both over- 
testing and undertesting. Nevertheless, it is vitally important to any 
vehicle-development program that, despite this fact, the most realistic test 
conditions available under the circumstances be sought, and the best tradeoff 
of factors influencing the degree of undertest ing or overtesting be selected. 
Arbitrary decisions concerning test techniques, levels, and durations can 
cause a considerable and sometimes an enormous unnecessary expense to the 
program in terms of reliability, schedule, and cost. For these reasons, little 
reliance should be placed on general specifications (ref. 146) purported to 
be applicable to every space vehicle. 

5.1.3 Acoustic Tests 

In the past, laboratory acoustic tests were used on an occasional basis to 
supplement qualification vibration testing of equipment. However, since two 
of the major causes of structural vibration are acoustic noise at liftoff and 
aerodynamic noise during the transonic and < l max regimes , it was a natural 

development to consider acoustic noise as a laboratory source of space-vehicle 
vibration. The Titan program was the first to utilize acoustic testing for 
the design development of vehicle structure subjected to acoustic noise 
(ref. 147). More recently, acoustic testing was used on the Apollo program 

(a) for qualification of equipment in lieu of qualification vibration testing, 

(b) for qualification of structure, and (c) for revision of vibration design 
and test requirements (refs. 148 through 152). Since the major objective of 
acoustic tests is to produce vibration in a large test specimen, it is 
common to call them vibroacoustic tests. 

5. 1.3.1 Test Facilities 

The facilities used for these tests are many and varied. For example, 
figure 36 shows the 0G0 spacecraft on the flatbed of a truck located in an 
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Figure 36 -Acoustic test of the OGO spacecraft near the discharge nozzle of a large blowdown wind tunnel. 


open field near the discharge nozzle of a large blowdown wind tunnel (ref. 153). 
Here, the acoustic level and spectra are dependent upon (a) the mass and 
velocity of the discharge flow, (b) the distance from the nozzle exit to the 
test specimen, (c) the angle between the exhaust direction and an imaginary 
line connecting the nozzle exit with the specimen, and (d) climatic conditions. 

A similar arrangement may also be used with a rocket engine as the noise 
source. However, if the acoustic test is intended to be performed as a 
‘‘free ride 55 on a wind-tunnel or rocket-engine test, a great discrepancy may 
exist between the desired and the actual test duration. In figure 37, the 
S-II stage of the Saturn V launch vehicle is shown during its installation in 
a large reverberant test facility (ref. 154). The progressive test facility 
used for the acoustic tests of the Apollo lunar module and the command and 
service modules is shown in figure 38 (refs. 150, 155, and 156). In the 
latter case, the external structure is excited by acoustic noise which is 
propagated down 16 ducts, each duct having its own noise source. The 16 sources 
may be operated and controlled together, independently, or in any combination. 
These noise sources, and those uaed in the previously mentioned S-II tests, are 
air modulators. These devices work on the principle of exhausting high-pressure 
air through an orifice whose cross-sectional area is modulated by an electro- 
magnet, which is controlled from an external electrical signal, usually a 
random-noise generator. In general, air modulators have limited spectrum 
range and control, usually from 50 Hz to 1 kHz, with a spectral maximum in the 
vicinity of 100 Hz. Noise generation above 1 kHz is usually determined by the 
4 ‘unmodulated 5 r flow noise cf the air through the orifice. Air modulators are 
commercially available in acoustic-power capacities up to 200 kW. 
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Figure 37 r Laboratory acoustic test of the thrust structure, aft skirt, and interstage of the S-ll 
stage, Saturn V launch vehicle in a reverberant test facility . 


5. 1.3*2 Comparison With Vibration Tests 

Although acoustic testing of space vehicles and large vehicle sections is 
still in its early development, a great many advantages have already been 
demonstrated over the vibration testing discussed in Section 5.1.2. If 
acoustic testing is performed on a large-enough section of the vehicle, so 
that the modal characteristics of the test specimen are similar to the modal 
characteristics of the flight configuration at the lowest frequency of interest, 
many of the following similarities with the flight configuration can often be 
achieved : 

1 . Similar (although not identical) vibration spectra at various 
locations throughout the specimen, except near where the 
specimen structurally interfaces with the acoustic facility. 

2. Antiresonant behavior is allowed to occur at the attachment 
of heavy structural members, including equipment items, when 
these members exhibit resonant behavior (refs. 138 and 139). 

3. Multidirectional resonances are allowed to occur. 

4. The effects on the integrity of the structure and the 
performance of the equipment from vibration occurring 
simultaneously in the three orthogonal directions can be 
determined. 
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Figure 38.- Laboratory acoustic test of the Apollo spacecraft in a progressive test facility 





5. 


The specimen is not restrained to an arbitrary amplitude or 
phase at locations where the vibration forces are applied. 

6. Combined vibration and acoustic loading can be applied to 
internal structures and equipment, which, under certain 
circumstances, can be made similar to flight conditions. 

7. If test conservatism is desired, it can be applied simultaneously 
in all directions and more uniformly over the frequency range. 

The major reasons for these improvements over vibration testing are (a) the 
ability to test a large and heavy structural specimen with a noise source of 
reasonable capacity, thus permitting the modal characteristics of the test 
specimen to be similar to those of the flight configuration over a wide 
frequency range; (b) the application of distributed loading over the external 
surface, rather than point-loading at one or more structural interfaces; and 
(c) fewer problems are encountered in selecting a reasonable acoustic-test 
spectrum than is found in selecting a reasonable vibration- test spectrum. 

5. 1.3.3 Test Requirements and Their Selection 


In spite of the bright picture just described, there are several important 
limitations in utilizing acoustic testing for demonstrating structural ade- 
quacy and equipment performance. At present, one of the two major technical 
limitations is the inability to select the proper acoustic test spectrum 
which will accurately produce vibration spectra throughout the test specimen 
which are reasonably similar to space vehicle vibration spectra at equivalent 
locations, before flight vibration measurements are made. This limitation 
applies to all three forms of acoustic testing presently being performed: 
free field (fig. 36), reverberant (fig. 37), and progressive wave (fig. 38). 
Equations (2a) and (2b) from Section 4. 1.2. 2 can be used to describe this 
problem. At liftoff, the vehicle structure is subjected to the rocket-engine 
acoustic noise that is propagated over the vehicle, as described in Section 3. 
If the specimen were tested in a reverberant facility, a joint acceptance 


2 

j^k would probably be found that differed from that at liftoff, 
to the difference in the cross-correlation functions (£,£*, t) 


owing primarily 
between 


liftoff and the acoustic test. 


A progressive-wave or free-field test might provide sufficiently close 
longitudinal correlation, but probably would have some significant difference 
in the circumferential correlation. In addition, most progressive-wave ducts 
develop transverse standing waves that can couple with the vibration of the 
external structure above a certain frequency. This limitation is even more 
restrictive when the acoustic-test spectrum is to be selected to provide the 
same vibration spectra as the vehicle structure exposed to aerodynamic noise. 
In this case, the propagation speed for the acoustic test will generally be 
different from the aerodynamic-propagation speed, which, in addition, is 
increasing during the launch -and -as cent phase. 
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Near the locations where the vehicle cross-sectional area changes, local high, 
fluctuating pressures from the separated flow during flight often cause a 
large gradient of fluctuating pressure over the external structure. However, 
it is unnecessary in most cases to provide the same gradient or the highest 
local fluctuating-pressure level during the acoustic test since equation (2a) 
shows that the pressure distribution is integrated over the area of interest. 
Besides, at low frequencies, the cross-correlation coefficient for aerodyanamic 
noise is much less than the coefficient for acoustic noise, thus producing the 
lower vibration shown in figure 20 at < 1. 

At high frequencies, and in the frequency range near the ring frequency of 
the external structure, coincidence occurs (i.e., the propagation speed is 
equal to the wave speed of flexural waves in the structure) , which causes an 
increase in the vibration response. Since the propagation speeds are 
different between flight and the acoustic test, the frequency ranges for this 
effect will not be identical (ref. 157). 

Another factor influencing the proper selection of the acoustic test spectrum 
is the difference in the acoustic radiation (i.e., air damping) between the 
vehicle during flight and the specimen during the acoustic test. Radiation 
is one of the forms of damping influencing the value of in equation (2b) . 

First, consider the space vehicle at liftoff. The external structure radiates 
sound externally to the free field (usually with little interference from the 
umbilical tower) and internally to the interior spaces. The external structure 
may be a payload shroud; a cylindrical or conical section enclosing a wide 
variety of equipment, as well as some primary and/or secondary structure; or 
a motor casing or tank for a solid or liquid propellant. Neglecting the 
latter, the external radiation from cylindrical shells is considerable for 
certain modes in the vicinity of the ring frequency, and for all modes above 
the coincidence frequency. References 53 through 56 describe the external 
radiation of unstiffened plates, cylindrical shells, and other surfaces. 

The internal acoustic radiation (for payload and equipment sections) is 
frequency -dependent (refs. 6 and 158 through 160). 

At low frequencies, the internal space usually acts as an air c ‘spring* 5 or 
stiffness, so that the internal radiation is negligible unless there is 
leakage or venting through small openings. At frequencies corresponding to 
the lower cc air modes’ 5 of the internal space, the internal radiation is 
usually small. Also, the air-mode natural frequencies and spatial distribu- 
tions are greatly influenced by large internal objects, such as the payload, 
equipment items, and other internal structures. In addition, these items 
also vibrate, contributing to the internal noise. At frequencies sufficiently 
above these lower air modes, the internal space is reverberant, and the 
radiation is dependent on the average absorption coefficient of the walls of 
the space and the objects. If the internal space and configuration of the 
test specimen is identical or similar to that of the flight vehicle, the 
internal radiation should also be similar (refs. 81 and 161). However, the 
external radiation of both may or may not be similar, depending on the type 
of test facility used. 
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If acoustic testing is performed in a free field, such as illustrated in 
figure 36, the external radiation of the test specimen and the flight vehicle 
should be identical. If acoustic testing is performed in a sufficiently 
large reverberant facility, so that the fundamental modes of the external 
structure are sufficiently above the lower air modes of the test facility, 
the radiation will be dependent upon the average absorption coefficient of 
the facility walls and the external surface of the test specimen, and thus 
will be considerably different from that of free field (ref. 160). If 
acoustic testing is performed in a progressive-wave facility, the radiation 
difference will be dependent upon frequency and the cross-sectional area of 
the duct (ref. 162). 

Second, consider the space vehicle in the transonic or regime and its 

influence on the acoustic radiation. The external structure radiates sound 
externally through the boundary layer into a partial vacuum, so that the 
amount of external radiation is reduced compared to that at liftoff. If the 
internal space is partially or completely vented, the internal radiation is 
also reduced. 

Another factor influencing the proper selection of the acoustic- test spectrum 
is the location of the microphone (s) for controlling and/or monitoring the 
test. Microphones placed near the test specimen may measure the radiation, 
of the external structure, as well as the intended loading or reference 
pressure spectrum. 

3. 1.3. 4 Assessment 

Although there are many complications in determining the effects of the 
differences in the loading and radiation damping, and of the microphone 
locations, it is possible to determine these effects and to calculate the 
necessary adjustment to the acoustic- test spectrum. However, this has yet 
to be done in practice. Also, it is anticipated that this adjustment would 
produce a desired acoustic- test spectrum that would be highly shaped. 
Presently available noise sources provide a relatively smooth test spectrum, 
though, so that some conservatism would probably be produced. In the future, 
narrow-band noise sources, or even shakers attached directly to the external 
structure, may be used to reduce the conservatism. 

Only the more obvious of the above factors are considered in the adjustment 
of the acoustic- test spectrum in testing performed to date. For example, a 
typical test may adjust for the differences in the joint acceptances simply 
by using the spatial-average flight-pressure spectrum for the acoustic test. 
Other factors are * ‘ considered 5 ’ by increasing the acoustic test spectrum 
until sufficient conservatism is achieved. This practice notwithstanding, 
the resulting conservatism in the vibration spectra, as measured throughout 
the test specimen, is usually much less than the conservatism obtained dur- 
ing conventional vibration testing. Thus , acoustic testing of fairly large 
sections of space vehicles has the reputation of providing more realistic 
vibration- test conditions without the gross overtesting normally associated 
with vibration testing. It is expected that in the near future the substitu- 
tion of acoustic testing for conventional vibration testing will be common- 
place throughout the industry. 
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Acoustic testing can play an important role in vibration prediction for 
equipment. This is especially important in the mid- and higher- frequency 
range, where (a) classical analysis cannot be used, (b) statistical-energy 
analysis may provide too much spatial or spectral averaging, and (c) extrapo- 
lation methods are too inaccurate. Figure 39 shows a typical comparison of 
flight-vibration spectra with the original vibration design and test spectrum, 
and a vibration spectrum measured at the same location during an acoustic 
test. If acoustic testing can be performed early enough in the program, and 
if dummy equipment is used, greatly improved vibration criteria may be pro- 
vided to equipment designers and subcontractors before the design is completed 
or qualification begins. 

There are some important management problems that should be considered. The 
initial cost of a large reverberant or progressive-wave facility, such as 
shown in figure 37 or 38, is extremely high. The operational cost should 
however be nominal. Thus, a large capital outlay will be required if these 
types of facilities are not available. If the noise source is a blowdown 
wind tunnel or a rocket engine, the operational cost of these facilities may 
be considerable unless the acoustic test is a ‘‘free ride’* on a wind-tunnel 
or engine test. In this latter case, the scheduling of the test may affect 
the acoustic- test schedule. Provisions must also be made for protecting the 
data-acquisition system from noise and climatic conditions. If operational 
equipment is included in the test, similar protection must be provided for 
the performance-monitoring equipment. 



FREQUENCY (Hz) 

Figure 39.- Comparison of vibration spectra from acoustic test with flight and with original 
design and test requirement. 
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If acoustic testing is used in lieu of vibration testing for the qualification 
of equipment, certain subcontract management problems could arise that probably 
would not occur if conventional vibration qualification testing were performed 
completely under the control of the individual subcontractors. It is assumed 
that the acoustic testing would be performed by the prime contractor, or by 
others under his direction, because the prime contractor would probably be 
responsible for furnishing the structural-test specimen, arranging for or 
furnishing the use of the acoustic- test facility, and managing the scheduling 
of the equipment items to be supplied by the several subcontractors. 

Problems of responsibility designation could arise if an equipment failure 
should occur during the acoustic test, and the subcontractor believes that 
the equipment item has been overtested or improperly monitored. Thus, 
sufficient operational, vibration, and internal acoustic measurements must be 
made so that equipment performance and the test environment may be adequately 
compared with the design and test criteria for the item. Also, the criteria ■ 
should provide for the application of the combined vibration and acoustic 
environments of the item. Other potential problems, such as the pretest 
adjustment of the test spectrum to which the equipment item is exposed, are 
less significant for acoustic testing than for conventional vibration testing. 
If it takes several test runs to complete the vibroacoustic qualification of 
the several equipment items, owing to item scheduling problems or the neces- 
sity for redesign and retest, fatigue failure may develop in the structural- 
test specimen because of the long exposure time. This result should not imply 
a structural inadequacy of the vehicle. Usually, minor repairs can be made 
to the specimen between test runs. Unless a major structural element becomes 
fatigued, the repair should produce only minor changes to vibration spectra 
throughout the specimen. 

Experience from the Apollo program has demonstrated that all these problems 
are nominal in comparison with the problems associated with conventional 
vibration testing, especially with its inherent overtesting to avoid the 
possibility of any undertesting. In addition, most major space vehicle 
contractors have large versatile acoustic facilities available for their use 
or have them in the planning stages (ref. 163). 

Acoustic or vibration tests are seldom performed on structures containing 
propellants, since solid-rocket motor cases arc usually unaffected by an 
external acoustic field, and liquid-filled tanks usually exhibit low vibra- 
tion from the mass- loading effects of the liquid at the lower predominant 
frequencies (ref. 5). 

5.1.4 Combined Environmental Testing 

As described in Section 4. 1.2. 3, failures may occur under combined environ- 
ments that may not occur when these same environments are applied sequentially. 
Since certain flight environments occur simultaneously, combined environmental 
testing is sometimes used for the laboratory qualification of space vehicle 
hardware. Various combinations of thermal, altitude (vacuum), sustained 
acceleration, vibration, and acoustic noise testing are most commonly used. 

In general, combined environmental effects are more difficult to predict on 
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equipment and its performance than on structure. For this reason, more 
combined environmental testing is performed on spacecraft and launch vehicle 
equipment than on launch vehicle structure. 

The additional cost of combined environmental testing is usually rather 
nominal, with the possible exception of sustained acceleration testing on a 
centrifuge. Centrifuge testing often requires rather elaborate test facilities 
and instrumentation, especially when combined with other environments. How- 
ever, these costs can often be justified when combined effects cannot be 
predicted on the basis of analysis and/or of sequential testing. The centri- 
fuge facility shown in figure 40 is the launch phase simulator. In addition 
to applying a combination of sustained acceleration, vibration, acoustic 
noise, and vacuum to small- and medium- size hardware, this facility can be 
programmed to simulate the time-dependent changes in these environments and 
thereby reproduce the proper environmental combinations usually found during 
the launch -and -as cent phase (ref. 164). 

5.2 FIELD TESTING 

In the field-testing category are various nonlaboratory tests, such as rocket 
engine and stage or payload static firings used to demonstrate the operational 
performance of various space-vehicle subsystems before flight. If good 
planning is done early in the program, these tests can also be used to pro- 
vide valuable data for predicting structural vibration or for revising 
earlier predictions. 

Solid- and liquid-propelled rocket engines are usually designed and tested 
quite early in the vehicle-development program, sometimes as a lead item 
even before the prime vehicle contractors are selected. Because of this 
early scheduling, these tests can be invaluable in providing acoustic- and 
engine-vibration data which can be used for early vibration predictions. 

This is especially important when a rocket engine incorporates new design 
features that can affect acoustic noise and vibration generation, but has 
not been tested or instrumented previously. This normally includes engines 
with new types of combustion chambers or thrust vectoring, or with different 
propellants or expansion ratios. However, the use of the rocket-engine 
acoustic-noise data must take into consideration any differences in the 
test- stand and launch-pad configurations that may affect the noise generation 
or transmission, such as the design of the flame deflector. 

Many of the components attached to the engine during the early firings are 
usually not of flight weight or configuration. In these cases, the spectral 
density G"(f) of the vibration applied to a flight component can be predicted 
from the spectral density G T (f) measured at the interface of the test-stand 
component with the engine, and from equation (26), where (f) , Z^(f), and 

Z^(f) are the mechanical (driving point) impedance of the engine, the 

test- stand component, and the flight component , respectively, as measured 
and/or predicted at the engine-component interface (ref. 105). However, most 
test- stand and flight components are sufficiently lightweight and small, so 
that test- stand vibration data obtained for these locations can usually be 
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Figure 40. - Laboratory test facility for combined sustained acceleration-vibration-acoustic-noise-altitude 
environmental testing of space-vehicle hardware. 



used for the flight components directly. Most rocket engines (and turbopumps 
for liquid engines) exhibit very high frequency vibration, so that it is not 
uncommon to measure vibration up to 20 kHz for these engines. The engines 
also transmit vibration mechanically to nearby sections of the vehicle. If 
they are included early enough in the test planning, vibration measurements 
can be made to determine the mechanical-vibration transmission from the 
engine(s). Mechanical -impedance measurements can be made and, together with 
the vibration measurements, used in equation (26) to predict the vibration 
of the space vehicle from this source. 

Although they usually occur late in the program, static firings of the various 
stages of the vehicle can provide valuable preflight data. Generally, only 
stages with liquid-powered rocket engines are statically fired. If the stage 
is the first stage of the vehicle, and the flame deflector and tie-down 
configurations of the test stand are the same as those of the launch pad, the 
acoustic noise and structural vibration for the static firing should be 
almost identical to the noise and vibration at liftoff, except perhaps at the 
top of the stage (where the interface with the next stage or payload will 
probably not be properly simulated). Thus, except for the top, the static 
firing should provide the best definition of lif toff ■ vibration throughout the 
stage. 

If the flame deflector and tie-down configurations of the test stand differ 
from those of the launch pad, the static- firing noise and vibration must be 
extrapolated to account for the difference. For locations far from the 
engine (s), where the mechanically transmitted vibration is negligible, the 
vibration prediction may be accomplished by one of two methods. First, for a 
detailed correction of the data, the vibration measured during the static 
firing can be scaled in accordance with the differences in the joint accept- 
ances and damping ratios, which must be calculated as described in Section 
4. 1.2. 2. Even the structural vibration from aerodynamic noise can be pre- 
dicted with reasonable accuracy when scaled in accordance with these differ- 
ences. Second, for a gross correction of the data, the vibration measured 
during the static firing can be scaled in accordance with the differences in 
the spatial-average fluctuating pressures for liftoff, or, with an additional 
correction such as that shown in figure 20 for vibration caused by aerodynamic 
noise. 

If the stage is an upper stage of the vehicle, the acoustic noise and struc- 
tural vibration for the static firing will seldom be the same as (a) that 
produced by the first-stage engines during liftoff, (b) that observed during 
the transonic and q^^ periods, or (c) that occurring during the powered 

flight of the upper stage, when in most cases, the only vibration is mechani- 
cally transmitted from the engines. At locations far from the engines, 
however, where the mechanical transmission is negligible, the vibration 
measured during the static firing can be scaled by methods described in the 
preceding paragraph to predict vibration of the upper-stage structure at 
liftoff, transonic, and q max periods. 

For locations near the upper-stage engines, it will be difficult to separate 
the relative contribution of acoustically- induced and mechanically- transmitted 


92 


vibration for the static firing. If knowledge of this relative contribution 
is required, special static-firing test facilities can be used if they are 
included in the planning early in the vehicle -development program. The 
following two test configurations must be implemented in sequence: 

1. The engine (s) must be isolated from the remainder of the stage 
by taking the engine- thrust loads out through the test -stand 
structure, which is separated and isolated from the test-stand 
structure that supports the remainder of the stage, while all 
other connections between the stage and the engine (s) are made 
flexible. 

2. The engine exhaust must be ducted away from the test stand for a 
considerable distance, so that the acoustic noise reaching the 
vehicle is negligible because of the great distance from the end 
of the duct to the vehicle. 

The first test configuration will permit the determination of acoustically- 
induced vibration from the upper- stage engine (s), with only negligible 
contribution from the mechanically-transmitted vibration, whereas the second 
configuration will permit the opposite. For locations near the upper- stage 
engine (s), the first configuration will allow the acoustic and vibration data 
to be extrapolated to liftoff, transonic, and q max conditions which occur 

during first-stage powered flight. The second configuration will allow the 
determination of mechani cally- transmitted vibration for upper-stage powered 
flight . 

It is not uncommon to use stage static firings as a form of vibroacoustic 
acceptance testing. The advantage of such testing is the performance evalua- 
tion of entire subsystems as they are exposed to combined structural vibration 
and acoustic noise, with the vibration occurring in three orthogonal 
directions simultaneously. The disadvantages are (a) in most cases, the 
acceptance vibration levels and durations are seldom found to be optimum, 
relative to the tradeoff between flight-adequacy demonstration and preflight 
degradation, as discussed in Section 5.1.2; and (b) major effect on the 
vehicle schedule may result if acceptance- test failures occur on the test 
stand, rather than in the laboratory. These advantages and disadvantages 
should be weighed carefully before static firings are designated for vibration- 
and acoustic-acceptance testing. 

5.3 FLIGHT TESTING 

The flight test is almost always used as the final demonstration of the 
adequacy of acoustic and vibration requirements. It is also the final 
demonstration of the structural adequacy and operational performance of the 
space vehicle. In case flight failure occurs, sufficient instrumentation 
must be onboard to determine the location of the failure and the probable 
cause. In too many cases, not enough instrumentation is provided, so that 
the failure location or probable cause is often not ascertained; then redesign 
must proceed on several fronts, usually at considerable cost and delay. Even 
with redesign, the real cause may not have been eliminated, which results in 


more flight failures, higher costs and delays, and sometimes contract 
cancellation or lost national prestige. In light of the history of flight- test 
failures, it is surprising to note that some space-vehicle development programs 
are still managed on a * ‘success -oriented 5 * basis (i.e., few contingencies 
allowed for flight failures). This philosophy usually changes, though, with 
the first flight failure. In many cases, there is not sufficient instrumentation 
because too little space or weight has been allocated for it. This problem 
can be avoided by program-management insistence that all instrumentation 
requirements be submitted during the preliminary design phase, when instrumen- 
tation space and weight are usually designated. Other problems arise when 
the telemetry subsystem, generally used in the monitoring of subsystem perform- 
ance, is inadequate. 

A large number of vibration and aeroacoustic measurements are usually needed 
during flight test to (a) demonstrate the adequacy of the vibration and 
acoustic requirements, (b) help determine the probable cause in case of 
flight failure, and (c) provide data for the design of future space vehicles. 
Measurement requirements usually far exceed the capacity of the telemetry 
subsystem - a situation that appears to grow worse with each new program. The 
main reason is that most vibration and aeroacoustic measurements require a 
wide bandwidth (usually from 20 Hz to 2 kHz, and sometimes up to 10 kHz) and 
a large dynamic range (usually from 30 to 70 dB) (ref. 165). Often more than 
a hundred temperature, pressure, dc voltage, strain, flow, low-frequency 
acceleration, and guidance functions can be telemetered over the same rf 
bandwidth used to transmit ‘one vibration or aeroacoustic measurement. Thus, 
the measurement requirements of the structural dynamicist are not usually 
received with enthusiasm by the instrumentation engineer. 

The situation has improved slightly by recent changes in the Inter-Range 
Instrumentation Group (IRIG) Telemetry Standards (ref. 166), which allows 
more high-frequency channels through constant bandwidth fm/fm and the exten- 
sion of proportional bandwidth fm/fm to higher frequencies. According to 
IRIG, single sideband am/fm and double sideband am/fm will be standardized in 
1970, thus providing a major improvement in telemetering dynamic measurements. 
However, a major breakthrough in telemetry technology will be required to 
provide a capability which will satisfy most flight vibration-measurement 
requirements . 

Other techniques can be employed to reduce the demand on the telemetry sub- 
system, including (a) onboard frequency analysis, where the rms or average 
output of a sweep or step-frequency analyzer is telemetered on a low-frequency 
channel; (b) onboard tape recording, where certain measurements are recorded 
during the periods of high vibration or aeroacoustic noise, and then 
telemetered later during a period of low vibration and noise; and (c) time 
division multiplexing, where certain measurements, usually from two to four, 
are sampled and telemetered in repeated sequence (e.g., 3 sec ‘‘on 5 ’ and 9 sec 
“off”). All these solutions, though, increase cost and onboard weight, and 
may permit some degradation in instrumentation reliability, so that all other 
possible solutions should be exhausted before these are considered. For 
example, in some cases a “boilerplate” spacecraft may be used during the 
early flight test of a launch vehicle. If the external aerodynamic 
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configuration of the boilerplate is identical to that of the actual spacecraft, 
aeroacoustic measurements should be made on the boilerplate and thus only 
vibration measurements will be required on the actual spacecraft. 

Two other major flight- instrumentation problems are those of calibration 
error and electrical noise, A calibration error can occur when the various 
gain and attenuator settings actually used throughout the data-acquisition 
and -reduction process differ from those scheduled or appearing in the 
instrumentation engineer’s or technician’s notebook. Thus, it is not uncommon 
to have one or more cases where, for example, 8-g rms are measured during one 
flight and 80-g rms during another flight at the same mission phase and at the 
same location on the vehicle - with the instrumentation group often providing 
6 ‘conclusive proof” that both measurements were accurate! 

Electrical noise in the form of “spikes’ 5 and random electrical voltages 
are also commonly observed on vibration and aeroacoustic measurements. This 
noise usually develops in the various signal conditioners, cables, or 
telemetry equipment, or in the ground recording facilities. Telemetry noise, 
is usually caused by poor rf signal strength with distance, or interference 
from rocket -engine -exhaust products or hypersonic plasma effects. Most 
structural dynamicists experienced with flight- test data can distinguish 
random vibration and aeroacoustic noise, from random electrical noise, except 
when the electrical noise is caused by an instrumentation component that is 
in itself sensitive to random vibration or aeroacoustic noise. Because of 
the high cost of flight vibration and aeroacoustic-noise measurements , it is 
certainly desirable to avoid these problems. For example, the likelihood of 
a calibration error, caused by an incorrect gain or attenuator setting, can 
be greatly reduced if an insert voltage or other “end-to-end’’ calibration 
is used, and if the ac calibration is applied immediately before flight and/or 
during periods of low flight vibration. The additional cost of “end-to-end’’ 
calibration equipment is usually negligible when compared to the overall 
cost of the measurements . Also, the “end-to-end’’ calibration signals may 
be used to set the gain and attenuation of the various instrumentation 
equipment, from the transducer signal conditioner in the vehicle to the 
frequency analyzer in the laboratory, thus reducing implementation costs. 

The main reason for acquiring flight-test data is the need for determining 
the adequacy of design and test requirements. However, before flight measure- 
ments are compared with these requirements, the flight data should be scaled 
to represent worst-case conditions. For example, at liftoff the vibration 
should be scaled according to the difference between the measured and the 
highest-rated engine thrust and the difference between the measured and the 
highest expected acoustic-noise transmission characteristics resulting from 
climatic conditions (wind, ambient pressure, temperature, and humidity). For 
transonic and Q max periods, the vibration should be scaled according to the 

difference between the measured and the highest expected density and aero- 
dynamic pressure. If sufficient measurements have been allocated to measure 
vibration at all locations of interest on the vehicle, some vibration and 
aeroacoustic measurements should be repeated on later flights to establish 
f light -to-f light variations. These should be considered when establishing 
worst-case conditions, unless the variance is so wide as to cause one to 
suspect the validity of the data. 


If the requirements exceed the flight data by a considerable margin, and if 
redesign of structure or equipment is contemplated for later vehicles (as is 
often done) , some weight reduction can often be realized by lowering the 
requirements without a significant sacrifice in flight adequacy or reliability. 
If the flight data exceed the requirements, then retesting and possibly some 
redesign should be performed to new requirements based on the flight data. 

If the flight data are to be used in predicting the vibration for a new 
vehicle, it will probably be desirable to divide both reference and new 
vehicles into zones and perform a statistical analysis of the flight data for 
each zone, as described in Section 4.3.2. Depending on the similarity between 
the reference and the new vehicles, one of the scaling methods of Section 
4.3.2, one of the frequency-response methods of Section 4.3.1, a modification 
of these, or a new method may be used for the vibration prediction. 

If the new vehicle consists of a new spacecraft on a standard launch vehicle, 
and if the vibration is transmitted to the spacecraft predominantly through 
the mechanical path (as discussed in Section 4.2.3), the spectral density 
G"(f) of the vibration applied to the new spacecraft can be predicted from 
the spectral density G* (f) measured at the interface of the launch vehicle 
with the reference spacecraft, and from equation (26), where Z^ (f) , Z^(f), 

and Z^(f) are t * ie mec ^ an i ca l (driving point) impedance of the launch vehicle, 

the reference spacecraft, and the new spacecraft, respectively (ref. 105). 
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OBJECTIVES OF SPACE-VEHICLE DESIGN DEVELOPMENT FOR VIBRATION 


Certain objectives must be achieved to ensure that the space vehicle is 
properly designed to withstand structural vibration occurring during the 
various phases of the mission, in sequence and/or in combination with other 
loads and environments, without structural failure or degradation. In addi- 
tion, design and test requirements for the various items of vehicle equipment 
must be selected to ensure their operational performance and structural 
integrity. To accomplish these objectives, the general tasks listed in the 
following paragraphs should be observed in every vehicle-development program: 

• Applied Loading. The characteristics of the fundamental sources of 
vibration and the vibration loading applied to the structure should 
be determined by analysis and/or test for periods of the flight 
when high vibration is expected. 

• Vibration Motion. The vibration motion of the structure due to the 
applied loading should be established by analysis and/or test. 

• Vibration Stress. The stresses resulting from the vibration motion 
should be determined for critical locations of the vehicle by 
analysis and/or test. 

• Allowable Stress. The allowable vibration stresses should be 
selected on the basis of the material properties, fabrication pro- 
cesses, time duration of the applied loading, and the effects of 
sequential and combined loads and environments (e.g. , static loads, 
mechanical shock loads, thermal loads, humidity, propellant atmos- 
phere) for critical locations of the structure. 

• Structural Integrity. The integrity of the structure should be 

evaluated by comparing the resulting stresses with the allowable 
stresses, including consideration of the following: (a) the statis- 

tical variation and accuracy of determining the applied loading, 
the vibration motion and stress, and the failure mechanism; (b) the 
consequences of accepting an inadequate structure or rejecting an 
adequate structure; (c) the reliability goal, if one has been set; 
and/or (d) the desired confidence in the integrity of the structure. 

• Equipment Requirements. Design and test requirements for vehicle 
equipment should be selected on the basis of (a) the determination 
of vibration motion due to the applied loading as outlined above, 
the determination of the equipment acoustic loading, the time 
duration of these loadings, and the accuracy of determining the 
applied loading, the vibration motion, and the equipment acoustic 
loading; (b) the statistical variation of the loading, motion, and 
the possible failure mechanisms; (c) the number of test specimens, 
the type of test to be performed (vibration or acoustic) , the degree 
of simulation, and the method of test control; (d) the consequences 
of accepting an inadequate equipment item or rejecting an adequate 


one; and (e) the reliability goal, if one has been set, and/or the 
desired confidence in the operational performance and structural 
integrity of the equipment. 

Documentation. All analyses and tests should be clearly documented 
and available to customer representatives for examination and 
review. Design and test requirements should be supported by 
documentation showing their derivation. 
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RECOMMENDATIONS FOR IMPLEMENTING VARIOUS DESIGN, ANALYSIS, AND TESTING 
TECHNIQUES 


All methods of design, analysis, and test described in Sections 4 and 5 have 
limitations of accuracy, cost, and/or timeliness to support program schedules 
and objectives. To provide a proper evaluation of structural integrity and an 
adequate selection of design and test requirements, in most programs it is 
necessary to apply several of these techniques during the various phases of 
vehicle development and to revise earlier predictions as later information 
becomes available. The selection of a particular technique is thus dependent 
upon the phase of vehicle development during which it is to be implemented. 

In this section certain techniques described in Sections 4 and 5 are suggested 
for implementation during specific phases of the vehicle-development program 
in order to satisfy the objectives listed in Section 6. 

7.1 PRELIMINARY DESIGN PHASE 

The evaluation of structural integrity should be initiated as soon as the pre- 
liminary vehicle configuration is selected; the initial determination of the 
structural characteristics is made on the basis of the static loading (engine 
thrust, vehicle acceleration, wind, etc.) and other loadings and environments 
(thermal, gusts, etc.). 

First, the applied vibration loading should be described. For liftoff, the 
acoustic loading from the rocket engine (s) can be determined from data sources 
referenced in Section 3.1. For vehicles with very 4 ‘clean 5 ’ external config- 
urations or for vehicle locations far removed from changes in the vehicle 
cross-sectional area, the aerodynamic fluctuating pressure loading at q max can 

be determined from data sources referenced in Section 3.3. To determine the 
applied loading from aerodynamic noise or buffet during the transonic period 
for any vehicle, and at q max for aerodynamically 4 ‘unclean 5 ’ vehicles, wind- 

tunnel tests must be used as discussed in Section 5.1.1. Usually, tests can 
be performed simultaneously with steady (i.e., static) aerodynamic-pressure 
measurements, so that close coordination between aerodynamicists and structural 
dynamicists should permit a considerable cost saving. Wind-tunnel tests are 
often performed during the early development phase of the program, but the 
detailed planning for these tests should be made during preliminary design. 

If rocket-engine development tests have been initiated, acoustic-noise measure- 
ments should be made to supplement or supersede the data referenced in 
Section 3.1. If not, planning should be made during preliminary design for 
the inclusion of acoustic measurements during the rocket engine test program. 
Also, vibration and mechanical impedance measurements to be used in the 
prediction of mechanically transmitted vibration should be made or planned 
during preliminary design, as discussed in Section 5.2. As an interim measure, 
vibration measurements from tests of an earlier engine should be used after 
they are scaled for differences in engine thrust, size, and other engine 
parameters. If the stage or payload is also subjected to static firing of 
its rocket engines, the applied loading for this condition should be 
determined. 
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Next, the vibration motion and stresses of the structure from the applied 
loading should be calculated. Since many or most of the structural details 
are lacking during preliminary design, it is often acceptable to represent 
the structure in these calculations by simplified mathematical models. This 
will permit an initial evaluation of the general adequacy of the overall 
structural design at little expense. This is especially important during 
early preliminary design when, in some programs, the structural design changes 
almost daily and a quick, inexpensive design evaluation is of great value. 

In the frequency range of the lower-order modes, the classical analysis dis- 
cussed in Section 4.1 should be used to perform the initial evaluation. 

However, simplified representations of the joint acceptance and the arbitrary 
selection of damping may be used to reduce calculation costs, and equivalent 
isotropic representation may be made of composite construction . Certain 
structural sections may be represented as flat or curved plates, cylindrical 
or conical sections, simple trusses, and lumped mass, or as their lumped- 
parameter equivalents. The contribution of the higher-order modes should 
be calculated by the statistical-energy analysis discussed in Section 4.2, 
but since the predominant stresses will probably occur in the lower-order modes, 
it may not be necessary to consider higher-mode contribution during preliminary 
design. Since structural integrity is dependent upon the effects of all com- 
bined and sequential loads on the structure, close coordination between loads 
engineers, stress analysts, and structural dynamicists is essential to a good 
initial design of the vehicle. 

Design and test requirements for vehicle equipment are usually needed as soon 
as the systems concept has been determined and preparations are made to solicit 
bids from subcontractors. If the requirements are to be submitted before 
the previously mentioned structural evaluation can be made, one of the extra- 
polation methods discussed in Section 4.3 and the test contingencies discussed 
in Sections 5.1.2 and 5.1.3 may be used to derive these requirements. As 
discussed in Sections 4.3.1 and 4.3.2, all extrapolation methods have serious 
deficiencies. Thus, it would appear advantageous to select more than one of 
these methods and use the most conservative prediction, since it is usually 
easier and less costly to reduce subcontractor requirements than to increase 
them. It is essential, however, that there be close coordination between the 
customer and the prime contractor’s structural dynamicists concerning the 
selection cf these requirements, so that the initial requirements are not 
‘‘cast in concrete’’ for the remainder of the program, but are modified as 
better information becomes available. If the design and test requirements are 
to be submitted during or after the initial structural evaluation, classical 
analysis and statistical-energy analysis may be used in the frequency ranges 
of the lower- and higher-order modes, respectively, to derive these require- 
ments in conjunction with the test contingencies discussed in Sections 5.1.2 
and 5.1.3. It is also suggested that little reliance be placed on ‘-general’’ 
specifications and other requirements that are not written or interpreted 
in terms of the specific mission or configuration of the space vehicle. 

Eesign-development tests to provide early information to the designer regarding 
the adequacy of his structural design should be planned during preliminary 
design. Close coordination between designers and structural dynamicists is 
required to ensure that the proposed tests are likely to provide the desired 
information. 
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During preliminary design, recommendations should be submitted to the customer 
for major tests to be implemented during later phases of the program. It 
is suggested that vibration and acoustic measurements be planned for static 
firings of stage or payload rocket engines, as discussed in Section 5.2, as 
well as for flight tests, as discussed in Section 5.3. For new rocket -engine 
configurations, consideration should be given to performing model engine tests 
to acquire external acoustic data to improve or supersede the initial static 
firing or liftoff acoustic predictions. Until the problems observed with 
vibration models are resolved, it is suggested that model tests not be used 
to predict vibration. 

Since the highest vibration usually occurs at liftoff and at transonic and 
Snax P er ^ oc ^ s » it i s highly recommended that acoustic tests of the vehicle or 

large vehicle, sections be performed to demonstrate structural integrity and 
to improve the vibration prediction, as discussed in Section 5.1.3. Since 
these tests are usually among the major ground tests of the program, close 
coordination between the customer and the prime contractor’s program manage- 
ment, structures, reliability, facilities, test, and structural dynamics groups 
is required, starting during preliminary design. As a major ground test, 
it is likely that it will be one of the main program constraints to flight 
testing. Also, it will probably be a major cost item to the program. 

If the vehicle or section structure is one of the first items off the produc- 
tion line, and good planning and test program control is implemented , the 
early acquisition of test results can be achieved. If failures develop, early 
redesign can be implemented. If no failures occur and sufficient strain 
measurements are made, and if vehicle weight reduction is desired, early 
redesign based on experimental data can be implemented. If vehicle equipment 
is installed, the acoustic test may also serve as a ccmbined vibration and 
acoustic-qualification test of the equipment items. If the test is performed 
before the actual equipment is available, dummy equipment should be installed. 
Vibration and acoustic measurements can then be used in the revision of the 
initial design and test requirements. In this case, it is suggested that two 
tests be performed: making use of dummies (1) to acquire early data for the 
revision of requirements, and (2) for the later qualification of the actual 
equipment . 

A ground vibration test, as discussed in Section 5.1.2, is not recommended as 
a substitute for an acoustic test because of the difficulty in determining a 
reasonable vibration test spectrum. However, if mechanically transmitted 
vibration from the rocket engine (s) is expected to be a problem, then a ground 
vibration test is recommended to determine the effects on the structure, with 
shaker (s) used as the simulated source (s) of vibration. Engine vibration and 
mechanical impedance measurements can be used with equation (26) to determine 
the test spectrum. The engine-vibration test progran should be planned during 
preliminary design. In certain cases, the same structural test specimen may 
be used for the acoustic test and the engine-vibration test. 
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It should be obvious from this discussion, then, that structural vibration 
planning should be started at the very beginning of the preliminary design 
phase. Also, it is apparent that sufficient flexibility should be allowed 
within the program to permit the application of improved definition of loads, 
vibration motion and stresses, failure mechanisms, and requirements to the 
vehicle design. 

7.2 EARLY DEVELOPMENT PHASE 

During the early development phase, great emphasis should be placed on the 
analytical effort. The evaluation of structural integrity, which was initiated 
during preliminary design, should be continued and refined during this phase. 
Since data from rocket-engine static firings, model-engine acoustic tests, and 
wind-tunnel tests provide better definition of the applied loading, and since 
the structural configuration becomes firm, the mathematical models and the 
classical analyses should be more detailed to provide improved vibration 
predictions. This is especially important if the earli r analyses showed 
that the vibration stresses approximated or exceeded the allowable stresses. 
Statistical-energy analyses of the structure should be performed and, together 
with the classical analysis, used to revise the initial design and test 
requirements for the equipment. Close coordination between designers and 
structural dynamieists will probably be required during this phase in order 
to resolve design problems as they arise. 

Design-development tests,' as discussed in Section 5.1.2, will probably be 
implemented during this phase. Structural dynamieists should be available to 
support the designer in the implementation of these tests and the interpreta- 
tion of the results. Assistance may also be required to evaluate acoustic 
tests and wind-tunnel tests in order to ensure that the data are adequate for 
subsequent use. Also, support may be required for the later implementation 
of stage or payload static firing and acoustic (or vibration) tests, such 
as resolution of problems concerning the use of facilities or the selection 
of instrumentation arid measurement locations. 

7.3 LATE DEVELOPMENT PHASE 

During late development, great emphasis should be placed on the experimental 
effort. Except for occasional revisions or refinements, the analytical evalua- 
tion of structural integrity should have been completed. Support to the 
designers concerning design problems and design-development test results should 
no longer be of major proportion. The major structural vibration tests are 
usually performed during this phase and will probably require a great deal of 
support from structural dynamieists . If proper planning has preceded them, 
ground acoustic (or vibration) tests and stage or payload static firings will 
probably provide a great deal of vibration and acoustic data that must be 
analyzed, validated, interpreted, and reported. Corrective action must be 
taken if failures occur. Test data should be scaled to account for any differ- 
ences between these tests and flight. Redesign should be made for any marginal 
structures. The design and test requirements for equipment should be revised 
based on static firing results and on new flight predictions. 
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During this phase, qualification tests of equipment are also performed. If 
feasible, it is suggested that most vehicle equipment be qualified for combined 
vibration and acoustic noise by installing the equipment in the structural test 
specimen for acoustic testing, as discussed in Section 5.1.3. If not, conven- 
tional vibration testing must be used for qualification. Support may then be 
required from the customer and the prime contractor’s structural dynamicists 
in the evaluation of test conditions, instrumentation, test results, and in 
the possible granting of waivers. 

7.4 FLIGHT-TEST PHASE 

Preceding flight testing, certain assistance may be required to verify to 
program management that with the successful completion of certain analyses 
and tests, the vehicle has demonstrated adequate structural integrity and 
that the equipment is qualified for vibration and aeroacoustic loading. Docu- 
mentation is often required to support the verification. With the completion 
of the preceding analyses and tests, this verification should be made with 
confidence . 

After each test flight, the flight vibration and aeroacoustic data should be 
reported and compared with the tests and analyses performed previously. If 
the mission parameters were not * ‘worst case’ 5 , the flight data should be 
scaled to worst-case conditions before this comparison. Structural integrity 
and the design and test requirements must be reevaluated, preferably before 
the next flight. Corrective action must be taken, as required, unless program 
management is willing to assume the risk of not incorporating these corrections. 

In case of flight failure, usually several possible causes are suspected. 
Vibration is often one of them. Adequate flight-vibration measurements should 
therefore be programmed to support possible postflight-failure analysis. 

Flight data should be carefully cataloged for later use, in case the space 
vehicle is selected for new missions or design changes are incorporated later 
in the program. 

Finally, papers should be written for publication in the technical journals 
to make the valuable data and experience acquired on the program readily 
available for use in future space-vehicle design. 
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CONCLUDING REMARKS 


The magnitude of the tasks outlined In the body of this report is expected to 
be considerable for most vehicle -development programs. For example, it is 
estimated that about 1.5% of the development ccsts for the Apollo spacecraft 
was associated with ensuring its structural integrity under aeroacoustic- and 
engine -induced vibration. Experience has also shown that failure to perform 
certain tasks can be a major source of risk to the successful development of 
the vehicle. Unfortunately, in many cases this risk can only be evaluated 
late in the program; that is, during flight test. 
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INTERNAL FORCE -AND -MOMENT VIBRATION- RE SPCNSE EQUATIONS FCR 
STRUCTURES EXPOSED TO AERO ACOUSTIC NOISE 

A. 1 FORCE-DISPLACEMENT AND MOMENT -DISPLACEMENT RELATIONSHIPS 

Most £ailure~prone space vehicle structures are those with structural section 
comprised of shallow shell elements (e.g., flat plates, cylindrical sections) 
such as shown in figure A-1. For these shell elements, reference 29 shows 
that the inplane force N^ , inplane shear force N^t bending moment , twist- 
ing moment and transverse shear force , (all per unit length) at 

location x in the X^Z plane is 


N 1 (x,t) 

“ S>[(u 1 + w/ap + v(v 2 w/a 9 )] 

(A- la) 

N 12 (x,£) 

= la>(i - v) (u 2 + vp 

(A- 1b) 

M 1 (x,t) 

= -K(w„ 1 + vw ) 

(A- 1c) 

M 12 (x,t) 

= D(1 - v)w 12 

(A- Id) 

Q 1 (x,t) 

= -D(W (11 Hr W 122 ) 

(A~1e) 


z 



Figure A-1. -Forces and moments on a shallow shell element. 
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where u, v, and w are the instantaneous displacements (at location x) in 

, X^, and Z directions, respectively; and a^ and a^, the radii of curvature 

in the X^Z , and X^Z planes, respectively. The partial spatial derivations 

are as follows: 


u^ = 8u/3x 1 
u 2 - 3u/9x 2 
v^ = 3v/8x 1 
v 2 ** 8 v/9x 2 


2 2 
w^ = 9 w/3x^ 

2 2 
w 22 = 3 w/3x 2 

? / 

W 12 ~ 9 w / 3x i 8x 2 

3 3 

w^ - 9 w/3x 1 

3 2 

W 122 = 9 w ^ 3x i 9x 2 


Similar equations may be derived for forces N 2 , and , and Q 2 , and for 

moments M 2 and in the X^Z plane. For an isotropic shell > the inplane 

and flexual rigidities are . 


% = Eh/ (1 - v 2 ) (A-2a) 

D = Eh 3 / 1 2 ( 1 - v 2 ) (A-2b) 


while for a honeycomb sandwich shell (ref. 167), with face sheets having the 
same material and thickness 


® “ 2E 2 h 2 (A-2c) 

D - 2E 2 h 2 h 2 /(1 - v 2 ) (A-2d) 


where E, h, and v are Young’s modulus, thickness, and Poisson’s ratio for the 
isotropic shell, respectively; E 2 , h 2 and v 2 , Young’s modulus, thickness, and 

Poisson’s ratio for each face sheet, respectively; and h^ , the half thickness 

of the core for a honeycomb shell. The instantaneous stresses are related 
to the various forces and moments by 


N 1 (x,t) 


+h* 

L 


S 1 1 (1 + z/a 2 >dz 


(A-3a) 
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+h* 


N 12 (x,t) = 


M-j (x,t) = 


M 12 (x,t) = 




L 


s 12 (1 + z/a 2 )dz 


h* 

+h* 


/ L 

* / _T_ -A. 


z(1 + z/a^)dz 


-h* 

+h* 


/ s ' 

^ T_ A 


•h* 

+h* 


/ -i 

-h* 


2 z(1 + z/a 2 )dz 


3 z(1 + z/a 2 )dz 


(A- 3b) 


(A-3c) 


(A-3d) 


(A-3e) 


For an isotropic shell, h* = 1/2; for a honeycomb shell, h* = h^ + h^. 
A. 2 FORCE-AND-MOMENT SPECTRAL DENSITIES BASED ON SPATIAL DERIVATIVES 


Similar to equation (A-2) the force-and-moment spectral densities for 
location x in the X 2 Z plane 

applied random loading are 


location x in the X 2 Z plane at each frequency f due to a spatially distributed 


G (x,f) =® 2 G <x,f) + v 2 G (x,f) 
1 U 1 V 2 


+ 23 


Jj^ 2 + (2v/a 1 a 2 ) + (v/a 2 ) 2 J G w 
( 

[ a i' + (v/a 2 ) ] + vG v, 




+ vG (x,f) 
u v 

12 [ 


G (x,f) = lg) 2 (1 - v) 2 |G (x,f) + G (x,f) + 2G (x,f) 

W 12 4 [_ u 2 V 1 U 2 V 1 


(A- 4a) 


(A-4b) 
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G M (x ’ f) 
1 

G M (x ’ f) 

12 

G n (x,f) 


d 2 r G 

L W 11 


(x,f) + v G (x,f) + 2 vG (x,f) 

w 22 w 11 w 22 


D 2 (1 - v) 2 G (x,f) 
W 12 


d 2 r G 

L w i 1 1 


(x,f) + G (x,f) + 2(3 (x,f.' 

122 111 122 


(A-4c) 


(A-4d) 


(A-4e) 


All spectral density terms on the right-hand side of equations (A-4) are 
expressible in the form 


N 


G t (x 5 f) = 


A 2 G (f) 

pr 


N 


2 2 

i=1 k=1 


<l> \ 00 ^ 00 H*<: f! )1 ^ (: f. ) j ^ 0 f ' ) 

(2 ”> 4f i4 M A 


(A- 5) 


where the relationships between the product of the mode shape derivatives 

[<j> /(x)4>^(x) ] and their spectral densities G (x,f) are given in table A- 1. 

1 ^ 

In this table, \J^(x),Y? ^(x) and (j)^(x) are the mode shapes in the i “ mode and 

in the and Z directions, respectively. Thus, with the appropriate 

substitutions of table A-I into equation (A-5) and the subsequent substitution 
of these into the equations of (A-4) , the spectral densities for the forces , 

and , and moments and ^ can be calculated. Spectral-density 

equations for forces ^ 2*\ 9 anc * ^2’ anc ^ moments M 21 t * ie P-^- ane 

can be similarly calculated. Rms forces and moments can then be calculated 
from equation (5b) , and rms stresses from the equations of (A-3) . 

It is apparent a considerable amount of calculation is required to determine 
forces and moments, which is unfortunate when the importance of these quanti- 
ties to the determination of structural integrity is considered. In certain 
cases, these calculations can be reduced somewhat by ignoring all terms in 

2 

equation (A-4) that are multiplied by v , since v ~0.3 for most aerospace 
metals. However, care must be exercised since other terms (i.e., not multi- 
2 

plied by v ) may be small at certain locations, and so cause appreciable 
un conservative errors. 

In most modal analyses, approximate methods are used to derive mode shapes. 
While an adequate approximate analysis can provide a reasonable estimate of 
the mode shape, a large error may be found in the spatial derivatives appearing 
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TABLE A- I. -RELATIONSHIPS BETWEEN SPECTRAL DENSITIES 
AND MODE SHAPE DERIVATIVES 


G t (x,f) 


^ (*)♦£(*) 


G (x,f) 
U 1 

G (x,f) 
2 

G (x,f) 
w 


[3^ i (x)/3x 1 ] [3i^ k (x)/3x 1 ] 


[3<^ ± (x)/3x 2 ] [3yj k (x)/3x 2 ] 

<j>i(x)0>k(x) 

[3iJj_ l Cx)/3x 1 ] <J> k (x) 


[3^ (x) /3x 2 ](j> k (x) 


[3i^ 1 (x)/3x 1 ] [3<,p k (x)/3x 2 ] 


[3iJ; i (x)/3x 2 ] [3^ k (x)/3x 2 ] 


[3^ (x)/3x 1 ] [3^ k (x)/3x 1 ] 


G (x,f) 

W 11 22 

G (x,f) 

12 


G (x,f) 

W 1 22 

G (x,f ) 

W 1 1 1 W 1 22 


[3^ i (x)/3x 2 ] [3 v?k (x)/3x 1 ] 
[3 2 4> ± (x)/3x*][a 2 .fr k (x)/3x*] 
[3 2 c|) i (x) /3x 2 ] [3 2 <J> k (x) /3x 2 ] 
[3 2 <|> i (x)/3x 2 ][3 2 <^ k (x)/3x 2 ] 
[3 2 (j) i (x)/3x 1 3x 2 ] [3 2 (fi k (x)/3x 1 3x 2 ] 
[3 3 ^> i (x)/3x 2 ] [3 3 <j> k (x)/3x 3 ] 
[3 3 (j> i (x)/3x 1 3x 2 ] [3 3 <J) k (x)/3x 1 3x 2 ] 
[3 3 1 1) i (x)/3x 3 ] [3 3 < t> k (x)/3x 1 3x 2 ] 
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in table A-I. Thus, in certain cases, forces and moments (and subsequent 
stresses) calculated by equations (A-4) and (A-5) may be in error by several 
orders of magnitude, especially for and , which are functions of 

second and third spatial derivatives (ref. 25). However, an alternate tech- 
nique, which is summarized in Section A. 3, can provide adequate estimates of 
vibration stresses even when approximate methods of modal analysis are used. 

A. 3 FORCE- AND-MOMENT SPECTRAL DENSITIES BASED ON EXTERNAL AND INERTIAL 
LOADING 

This technique is an extension of the modal displacement and modal accelera- 
tion methods of reference 24. The instantaneous internal force or moment P 
at location x due to instantaneous forces F applied at locations x is 

P(x,t) = ^ B p (x,x)F(x,t) (A- 6 a) 

where B p (x,x) can account for the effects of redundance if the structure is 
redundant (ref. 25). The force F has two components: 

F(x,t) = F e (?,t) - F..(x,t) (A- 6b) 

where F^ and F.. are the external force and the inertial force, respectively. 
For a distributed loading on a distributed structure: 

P(x,t) = Tb (x,i)p(i>t)dl 

A J 

+ J B (x,x)p (x)h(x)ft(x,t)dx (A- 7) 

where p(£,t) is the instantaneous external pressure at location g. For an 
isotropic shell of uniform thickness, the surface density is 

[p (x)h(x) ] = ph (A-8a) 

whereas for a sandwich shell 

[p (x)h(x) ] = 2(p 1 h 1 + P 2 h 2^ (A- 8b) 

Under random external loading, the force or moment spectral density is 

G p (x,f) = G p (x,f) + G p (x,f) + G p p (x,f) + G p p (x,f) (A-9a) 

P W P w w p 
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where 


G p (x,f) = 
P 

G (x,f) = 
w 



A 



5 f)d5dg_' 

Bp( x »E)Bp( x >2.’ ) tP (x)h(x) ] [p (x 1 )h(x’ ) ] 


(A- 9b) 


(2irf)^G w (x,x’ ,f)dxdx’ (A“9c) 

and where G_(g,?', f) is given in equation (3) and G (x,x'f) is given in 
r W 

equation (2), substituting [<|>^(x) (^(x* ) ] for [<}k (x^^Cx) ] • It is reasonable 

to assume that the external loading and the response in any mode is uncor- 
related, so that 


G p p (x,f) = G p p (x,f) = 0 (A-9d) 

p w w p 


Since P may be any force or moment, such as those expressed in equations 
(A-l) and (A-3) , the force and moment spectral densities may be obtained 
from repeated solutions of equations (A-9). For example, the bending -moment 
spectral density at location x in the X 9 Z plane can be found when B is 


substituted for B p : 


G M (x ’ f) 
1 


JJ B (x,pB (x,£')G a,r,f)dCdi' 

A 1 1 

+ JJ B (x,x)B (x,x T ) [p(x)h(x)][p(x l )h(x ! )] 
S 1 1 


(27rf) 2 G^(x,x ! ,f)dxdx* (A-l 0) 

The transfer function B p between the internal forces or moments, and the 

applied force, given in equation (A-6a) , can be obtained from static-stress 
analysis, and in fact is normally calculated by stress analysts for static- 
load analysis* However, the magnitude of the task is still considerable 
because B p must be obtained for every combination of potential failure loca- 
tion x and location x of the structure (or location £ of the structure exposed 
to the external-pressure field)* Thus, it is apparent that the calculation 
of stresses is considerably more complex and expensive than the calculation 
of displacements or accelerations. 
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A. 4 ASSESSMENT 

Examination of equations (A-9) shows that they contain no spatial derivatives. 
Thus, when approximate methods of modal analysis are used, appreciably greater 
accuracy may be expected from the use of equations (A-9) than from equations 
(A-4), especially for the calculation of the spectral densities for moments 
, M^, M^, and , and forces and Q^, which are functions of second and 

third derivatives. 
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